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ABSTRACT 


This report evaluates the performance of five digital range 
rate filters. A range rate filter receives an input of range data from 
a radar unit and produces an output of smoothed range data and its 
estimated derivative range rate. 

A comparative study of the filters is made through simulation 
on an IBM 370. Two of the filter designs are also implemented on a 
6800 microprocessor-based system. 

Comparisons are made on the bases of noise variance reduction 
ratios and convergence times of the filters in response to simulated range 
signals. 
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I . INTRODUCTION 


The range rate filters discussed in this report are digital fil- 
ters designed to accept a digitized input from navigational equipment, 
representing a range value, and to output a smoothed range value and its 
estimated derivative, range rate. Five filter designs are evaluated. 

Each is based on different design concepts; hence, the filters have 
widely varying capabilities. 

The purpose of this report is to present a comparative study of 
the five digital range rate filters. With this In mind, a section 
detailing the theory of operation for each filter is provided. As a 
comparative study, the main objective is to present actual test data in 
such a manner that direct comparisons can be easily drawn between the 
widely varying filter types. Based on this data, conclusions are made 
as to filter perforitiance in terms of variance reduction ratios, conver- 
gence times and ease of implementation. 

The report is divided into two major sections — full-scale 
computer simulations and microprocessor simulations. In the first section, 
the five filter designs are implemented on an IBM 370 computer. Given 
the higher precision available on a full-scale computer, these simula- 
tions represent best case results. Each filter is subjected to the same 
tests performed at two sampling rates: 2500 samples/second and 6.25 

samples/second.. The data is presented in tables for each filter. 

From the information obtained from the full-scale computer simu- 
lations, two filter designs are chosen to be implemented on a micropro- 
cessor based system. The 1 i. tions on speed and word size imposed by 
the microprocessor require special attention be given to programming the 
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filter algorithms. The resulting simulation data is presented along with 
equivalent IBM 370 simulation data for direct comparison. 



II. DESCRIPTION OF FILTERS 


A. Filter 1 

Filter 1 is a digitized version of the continuous system des- 
cribed by the transfer functions 


w - M.- sKi + K2 

- RM -■ + sKl + K2 


( 1 ) 


Hr^(s) 


ROOT _ sK2 

P + sKl + K2 


( 2 ) 


for range and range rate, respectively. The block diagram for this analog 
system is shown below. 



FIGORI 1 - Filter 1 (continuous system) 

The discrete system is constructed by replacing the integrators in Figure 
1 with discrete integrators using Euler's approximation. A delay must 
also be added to the feedback loop to make the system physicalTy reali- 
zable. This method of mapping from the s-pTane to the z-plane is 
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5 

discussed in further detail in Appendix I. The resulting discrete 
system is presented in Figure 2. This system is characterized by the 
transfer functions: 

.. _ RE _ KIDT + K2DT“ - z'MKIDT) 

M " 1 + 2“M-2^“KiE)T + K2DT^} V^-ni - IQDT) (3) 

.. _ RD0T _ K2QT - z'^K2BT) 

■ "W ■ 1 + z-^l-^ > KIDT + K2DT^} + z^^U - KIOT) (4) 


The classification of filters described by equations (1) and 
(2) is termed constant gains. This type of filter is simulated for five 
sets of parameter pairs, K1 and K2. In addition, if the parameter pairs 
are chosen such that 


K1 = b/Z 


(5) 


and 


K2 = b^ 


( 6 ) 


where b is a frequency scaling parameter, the constant gains filter 
becomes a bandlimited differentiator.^ Filter 1 is also simulated 
for five parameter sets chosen in accordanee with equations (5) and (6). 

B. Filter 2 

Filter 2 is an alpha-beta tracking filter described by the 
difference equations; 
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RE(n) = RDI(n) + a (RM(n) - RDI(n)) 


(7) 


RDOT(n) = RDOT(n - 1) + ^RM(n) - ROI(n)) (8) 


RDI(n + 1) = RE(n) + DT{RDOT(n)) (9) 

A block diagram is presented in Figure 3. 



FIGURE 3 - Filter 2 

At the sample period three values are computed ; the estimated 
range at that sample period, RE{n); the estiimated range rate at that 
sample period, RDOT(n); and a predicted value of the range at the 
(n + sample period, RDI(m-l). The input value is RM(n), and 
DT iS' the sample interval . The response of this type filter is 
con trolled by the parameters a and 8. For this appli cation, the re- 
lation 





(10) 


6 = 



is used to yield the maximufh maneuver- following capabilities for a 
given level of noise rejection. This results in the trade-off between 
noise rejection and maneuver- foil wing capability being controlled 
by only one parameter, a. 

The alpha-beta tracker remains a stable system for a and g . 

3 

bounded by 


a > 0, @ > 0 and 2a + $ < 4. 


( 11 ) 


For a = ^ = 1 ( a point on the curve denoted by (10) within the region 
specified by (11)), the system has a wide bandwidth and reacts quickly 
to transients. However, the transient-following capability is offset 
by poor noise rejection. The converse occurs for a « 1, where the 
system is slow and has a comparatively low bandwidth. Noise rejection 
is much improved, but the system reacts more slowly to transients. 

The transfer functions for the alipha-beta tracked are 


Hr(z) 

H|^j^(z) = - 


_ M _ Qt (l-a)z"^ 

RM 1 + (-2 + a + + (l - a)z"^ 

[DOT ^ 1 i - gz'^ 

RM “ DT 1 + (-2 + a + $)z“^+ (1 - a)z-^ 

RDI _ (a + S)z”' - az”^ 

RM " I ¥ (-2 # a + (1 - a)z^^ ’ 


( 12 ) 


(13) 


( 14 ) 
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where Hn is the transfer function for the predicted range, RDI, 

Kp 

It is interesting to note that the transfer function for Filter 1 
can be made identical to for the alpha - beta tracker if the para- 
meter pair K1 and K2 are chosen such that 


a = K1 * DT 


(15) 


and 


$ = K2 * Dt" , 


(16) 


the range transfer functions^ for Filter 1 and 


P 


for 


the alpha-beta tracker, are also very stmilar in that they differ 
only by a factor of The main difference between the two filters 
is that the transfer function coefficients are dependent on sampling 
rate in Filter 1 and independent of sampling rate in Filter 2, the 
alpha-beta tracker. Simulation verifies that the two filters can 


be made to perform identically in the range rate mode and very closely 
in the range mode. 


G. Filter 3 


A functional block diagram of Filter 3 is presented in Figure 4. 
Eighteen registers containing the most recent and the previous seven- 
teen range inputs are used in a cascaded simple average smoother to 
derive range rate. The sum of the six oldest range yaTues is sub- 
tracted from the six most recent values. The result is a sum of six 
terms in which each term is the Change in range that occurs over 
twelve sample periods. Division by six yields an average change in 













i 

i 

1 . 


I 

t. 


j / 
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range over twelve sample periods. Further division by 12 * DT then 
results in an estimated range rate value. 

An N sample smoother is used to determine smoothed range. The 
sixteen most recent range inputs are averaged to yield a smoothed 
range delayed by 7.5 sampling periods from the current sampling 
period. The current range rate is used to determine a correction 
term that is added to the smoothed range to update it to the current 
sampling period. 


D.. Filter 4 


Filter 4 is a digital tracking filter constructed by decreas- 
ing the number of storage registers and adding feedback to the design 
of Filter 3. A block diagram of Filter 4 is presented below. 



FISURE 5 - Filter 4 

In this design the values contained in the four storage registers 
now represent range error terms, rather than range terms, due to the 
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addition of feedback. An average of these four registers is taken 

resulting in an average range error term, this value is added to the 

previous estimated range to yield the new estimated range. The 

estimated range rate value is determined by dividing the average range 

error term by the sampling period DT. The estimated range is then 

fed back to the sampled range input to produce a new range error term. 

Four is the maximum number of storage registers (samples) 

4 

that will result in a stable system. 

E_. Filter 5 


Filter 5 is a tracking filter similar to Filter 4 except for 
the addition of skip sampling. In this design, although the sixteen 
storage registers are still updated every sample instance, the range 
feedback term is updated only every sixteen sample periods. Estimated 
range and range rate are also available only at sixteen^sampl e-period 
intervals. A block diagram of Filter 5 is provided in Figure 6. 

Every sixteenth sample^two range values are determined: the 
estimated range at that sample instant and the feedback range which 
lags the estimated range by 7.5 sample periods. The feedback range 
is subtracted from each of the next sixteen range inputs so that the 
value in the sixteenth storage register represents the change in range 
over 8.5 sample periods, and the value in the fifteenth refister re^ 
presents the change in range over 9.5 sample periods. This same pattern 
Gontinues through the. first storage register which contains the change 
in range over 23.5 sample periods. An average taken of the sixteen 
storage registers then yields a value that represents the smoothed 
change in range over sixteen sample periods. This value is added to 















the present feedback range to determine the feedback range to be used 
for the next set of sixteen range inputs. The smoothed change in 
range over sixteen sample periods is then multiplied by the ratio 
7.5/16, and the result is added to the new feedback range to yield 
the estimated range for the present sample. The estimated range rate 
can be derived by dividing the smoothed, change in range by sixteen to 
determine the change in range over one sample period. Division by the 
sample period then yields the estimated range rate. 



III. IBM 370 COMPUTER SIMULATIONS 


A, DesGription of SImuTatTnq Program 


Program FLT is the FGRTP.AN program used in the I3M 370 
simulations. In this program an initial range value and range rate 
value are entered to create a range signal of constant slapeo This 
ramp simulates an object moving with constant velocity and is termed 
the true or actual range. A random number generator is used to create 
an array of normally distributed noise values of specified mean and 
standard deviation. Each filter is simulated for three range inputs. 
Simulation RUN 1 tests the filter response to the ramp, RUN 2 tests 
with the noise and RUN 3 uses a range input that is the summation of 
the ramp and the noise. 

Filter calculations are performed by subroutine CALC. This 
subroutine also determines the input range value for a given sample 
period. That value is entered into the filter algorithm, and an 
estimated range and range rate value for that sample period are deter- 
minedo The value IDI in the program deteirmi nes how many times this 
process is repeated before program control returns to FLT. CALC returns 
to FLT the value at the most recent sampling period of the actual range, 
the range input, the estimated range, the actual range rate and the 
estimated range rate. These values along with the range error and 
range rate error, as determined from the differences of the estimated 
and actual range and range rate values, are output in tabulated form 
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and saved in an array to be used in plotting. The value IDRI 
determines how many times subroutine CALC is called. 

At the end of a simulation run, program FLT calculates • 
statistics on the steady state values of various range and range rate 
signals. The subroutine BECORI is called from the International 
Mathematical and Statistical Library for this function. The mean and 
standard deviation are calculated for a set of sample values of the 
input range, 'estimated range and the estimated range rate. For all 
simulation runs at 2500 samples/ second and 6,25 samples/second, the 
last 50 sample values of each of the three signals are saved in an 
array to be analyzed by BIGORI. Since the last 50 samples are used, 
it is necessary that the simulation runs be of sufficient duration 
to allow initial transients to abate. This will insure that the samples 
represent steady state signals so that statistically accurate results 
can be obtained. 

Program FLT contains two plotting routines for graphical 
representation of data. The first, PLOT B, implements the printer 
as a plotter. Although it provides limited resolution it is useful 
as a quick check of simulation results. The second routine, MCPLOT, 
formats data for a drum type plotter. After simulation RUN 3, four 
plots are generated for a particular filter - two range plots and 
two range rate plots . The first plot of each pair contains the 
graph of the actual signal and the resulting estimated signal • 
over the full duration of that simulation fun. however. In some cases 
the magnitude of the initial transient causes the scale of the plot 
to be too large to display any information on the steady state signal. 
For this reason, the second plot, an after-transient graph of the signal. 
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is made. It displays the actual and estimated signals after a time 
delay to allow for decay of the transient. These four plots are 
generated for every filter being simulated. At the end of program 
FLT, after all filters have been simulated, four additional plots are 
generated. These contain the transient and after- transient graphs of 
the range error and range rate error signals. As an example, if 
two filters have been simulated, each of the four final plots will 
contain the graph of two signals one for each filter. This allows 
for a comparison of filter results. 

As a clarification, only one type filter can be simulated 
during a given execution of program FIT, where type of filter refers 
to Filter 1, Filter 2 , etc. References made in the above text to 
different filters being simulated during the same execution of pro- 
gram FLT mean different parameter pairs such as K1 and K2 used 
in the same type filter. 

A listing of program FLT, subroutine CALC for each of the five 
filter types, and the plotting routines, PLOT B and MCPLOT, are given 
in Appendix A. 


B. Test Cjanditions 


As previously noted, three simulation runs are made to evaluate 
each filter. For the first run, the input signal is a ramp function 
representing actual range and is described by the equation: 


r(t] = 60800 - 28t. 


(17) 



This is a ramp of negative slope corresponding to an initial range of 
18.532 Km (60800 ft.) and a range rate of -8.5344 m/seCo(-28 ft. /sec.). 
The purpose of this run is to evaluate the ability of the filter to 
accurately track an object moving at constant velocity. This input 
signal has been used in previous investigations of range rate filters. 

The second simulation run determines how well the filter re- 
jects random noise. The input signal is Gaussian noise with a mean 
(p) of zero and a standard deviation (a) of 30.48 m (iOO ft.) 

The purpose of the third run is to evaluate the response of 
the filter to noisy range data. The input signal is a summation of 
the inputs for the previous two runs. This Input can be termed a 
noisy ramp. This run shows what effect the addition of noise has on 
the ability of the filter to track the actual range value. 

In the first set of tests, all filters are run at a sampling 
rate of 2500 samples/seeond for a period of 75 seconds. Results are 
taibulated at one second intervals, and the steady state statistics, 
mean and standard deviation, are taken on the last 50 samples. 

The second set of simulations are made at 6.25 samples/second 
for a duration of 192 seconds. Results are tabulated every 16 sample 
periods or 2.56 seconds. Steady state stai,isties for these simulations 
are taken on the last 50 samples. 

G,. Test Results 

Data for each f i 1 ter i s di vi ded i nto two secti ons : that ob- 
tained at a 2500 samples/seeond sampling rate and that obtained at 


6.25 samples/second. Data for each filter is provided in tables and 
in many cases also presented in graphical form to aid in visualizing 
actual range and range rate signals. A summary of the filter's per- 
fomance is given at the beginning of the data for that filter. 

Four tables are presented for each filter. The first table 
gives the input and output range standard deviations for each of the 
three simulation runs. The second table gives the output range rate 
standard deviations for each of the three simulation runs. These values 
give a measure of a signal's variance from its mean. In the case of 
simulation RUN 1, where the range input is the ramp representing 
actual range, the desired range output signal is an identical ramp. 
Hence, the standard deviation of the input and output range signals 
should ideally be equal. Since the slope of the ramp is constant, 
it is expected that the estimated range rate signal produced by the 
filter should also be constant and hence its standard deviation will 
approach zero. During simulation RUN 2, the input range signal is 
(Saussian noise of zero mean and known standard deviation. If the noise 
is indeed attenuated by the filter, then the estimated range output will 
show a reduction in standard deviation from that of the input range 
signal. Since the input represents a constant range value (namely 
zero range) with noise added, the optimum estimated range rate signal 
returned by the filter would again have a standard deviation approach- 
ing zero. The range input provided to the filter in simulation RUN 3 
is that of the Gaussian noise superimposed on the ramp. If both the 
ramp and the noise are independent, random signals and the filter is 
a linear system, then 



19 






+ N 


(18) 


where a| is the variance of the estimated range output for the ramp 
input, is the variance for the estimated range output for the noise 
input and ^ is the variance for the estimated range output for 
the ramp plus noise or noisy ramp input. In this case a|, 
and O’! + would represent the square of the estimated range standard 


deviations for 


1, RUN 2 and RUN 3, respectively. However, although 


the ramp and noise are independent, the ramp is not a random signal 
so the above relation does not hold true. For this reason, RUN 3 
is necessary to determine the results for a noisy ramp range input. 

The third table contains the maximum steady state range and 
range rate errors for the three simulation runs. This information is 
presented to give another measure of the range over which an estimated 
range or range rate value will deviate from the actual value. 

The fourth table lists the range variance reduction ratio. 


KR, and the range rate variance reduction ration, KRR. These 


ratios are determined frcmi the relations: 


i(D - variance of range output 
vari anee of range i hpt 


(19) 


and 


KRR = variance of range rate output 
“ variahee of range ihput 


(20) 
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The variance terms are the square of the standard deviations presented 
in the first and second tables under RUN 2= This is the simulation 
run conducted with Gaussian noise as the range input. Whereas the 
standard deviations provide information on specific signals, the 
reduction ratios offer a standardized method for comparison of filter 
results. 

The fourth table also lists convergence time x. This was 
arbitrarily defined as the time required for the estimated range to 
reach a steady state condition such that the magnitude of the range 
error does not exceed 3.048 m (10 ft.). The data used to determine 
T is taken from simulation RUN 1, where the range input is the ramp. 

1, Filter 1 

Simulation data at 2500 samples/second for Filter 1 
is presented in Tables I-IV for the constant gains filter and in 
Tables V-VIII for the bandlimited differentiator. Data at 6.25 
samples/secohd is presented in Tables IX -XI I for the constant gains 
filter and in Tables XIII-XVI for the bandlimited differentiator. 

A graphical display of the data at 2500 samples/second is shown 
for two sets of parameter pairs (Kl, K2) for the constant gains 
filter in Figures 7-12 and for the bandlimited differentiator 
in Figures 13-18, Data at 6.25 samples/second is shown for the 
constant gains filter in Figures 19-24 and for the bandlimited 
differentiator in Figures 25-30. 

The data indicates that Filter 1 has higher noise re- 
jection at the higher sampling rate. Since this filter is a 
digital approximation to a continuous system, this is expected. 
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The digital system more closely approximates a continuous system 
as the sampling rate increases. 

For the band! imi ted differentiator (BID) the effect of 
b, the frequency scaling parameter, can be noted by observing the 
variance reduction ratios. In the band! imi ted differentiator, the 
cutoff frequency is determined by the frequency scaling parameter. 

As the scaling parameter is reduced, the bandwidth of the system 
decreases. Hence, more lower frequency components of the noise 
are attenuated and the variance reduction ratios decrease. This, 
however, also reduces the system's responsiveness to higher frequency 
changes ‘in range signals which results in a reduction in maneuverr 
following capabilities. This is shown by the increase in conver- 
gence time as the frequency scaling parameter is reduced. 

A comparison of the variance reduction ratios shows a 
slight advantage for the bandlimited differentiator over the constant 
gains filter. While the constant gains filter produces good range 
and range rate data, there appears to be a trade off between the 
two. With the bandlimited differentiator, changing the frequency 
scaling parameter to improve range data also has the effect of 
improving range rate performance (at the expense of convergence 
time). 

There are a variety of methods to derive a discretized 
approximation to a eontlnuous system. The method used in Filter i 
produces a system that is a predictor, i.e., the estimated range 
returned at the n^^ sampling instant is a predicted value of the 
range at the (n + 1) sampling instant. This effect can be ob- 
served in the maximum steady state range error for simulation 



TABLE I 


FILTER 1 RANGE STANDARD DEVIATIONS(FT. ) AT 2500 SAMPLES/SEC. 



RON 1 

. input = ramp 

RUN 2 

input = noise 

RUN 

inout = rain 

3 

+ noise 

FILTEK 1 

constant gains 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STN0. 0EV. 

INPUT 

RANGE 

STNO. DEV. 

OUTPUT 
RANGE 
STND. DEV. 

INPUT 

RANGE 

STMD. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

K1 = 0.7071 
K2 = 0.4‘S 

G.1634 

0.1634 

110.7364 

0.8773 X lO"^ 

■ 

110.7339 

0.2463 

Ki = 0.75 
K2 = G.125 

■ 

0.1634 

0.1634 

■ 

110.7364 

■ 

0.9391 X 10’^ 

110.7339 

0.2523 

Kl = 0.7 
K2 = 0.3769 

0.1634 

i 

^ 1 

0,1634 

110.7364 

0.8 717 X 10'^' 

110.7339 

0.2459 

Kl = 0.49 
K2 = 0.98995 

: 

0.1634 

0.1632 

11G.7364 

0.5529 X 10"^ 

i 

110.7339 

G.2142 

Ki = 0.75 
K2 = 0.45 

G.1634 

0.1634 

110.7364 

0.9325 X lO"^' 

110.7339 

0.2511 


ro 

rc 































TABLE II 


FILTER 1 RANGE RATE STANM DEVIATIONS ( FT./SEC. ) AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


input = range 

Input = noise 

iinput = ramp + noise 

FILTER 1 

OUiTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 

eeNSTANT GAIiKS 

STND. DEV. 

sTnd. dev. 

STND. DEV. 

K1 = 0.7071 

9.0 

0.5652 X 10'^ 

0.5652 X 10"^ 

K2 = 0.4'5 




K1 = Q.75 

0.0 

0.1575 X 10"^ 

0.1575 X 110'^ 

K2 = 0.125 




K1 = 0.7 

0.0 

0.4740 X 10"^ 

0.4740 X 10‘^ 

K2 = 0.3769 




11 

o 

• 

0.0 

0.1238 

0.1238 

K2 = 0.98995 




Kl= 0.75 

0.0 

0.5662 X lO"^ 

0.5661 X 10“ 

K 2 = 0.45 




















TABLE III 


FILTER 1 mmm RANGE{FT.) and range RATE{FT./SEG.) errors at 2500 SAMPLES/SEC. 



RUN 1 

RUN 

2 

RUN 3 


input 

= ranip 

input = 

flQise 

input = nanip + noise 

FILTER 1 

MAX 

MAX 

MAX 

MAX 

MAX 

MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

CONSTANT GAINS 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

K1 = 0.7071 

0.112 X 10“^ 

0.14:0 X 10“^ 

0.278 

0.323 

0.290 

0.323 

K2 - 0,45 







K1 = 0.75 

0.108 X 10“^ 

0,221 X 10"^ 

0.223 

0.138 

0.234 

0.139 

K2 = 0.125 







K1 = 0.7 

0.112 X 1!0"^ 

0.225 X 10"^ . 

0.247 

0.264 

0.258 

0.264 

K2 = 0.3769 







K1 = 0.49 

0.106 X 10'^ 

0.464 X 10“^ 

0.322 

1.395 

G.332 

1.395 

K2 = 0.98995 






i- 

Ki = 0.75 

0.112 X 10“^ 

1 

0.241 X 10"^ 

0.250 

0.329 

0.261 

0.329 

K2 = 0,45 






j. 


ro 

4 =« 


TABLE IV 


FILTER 1 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 2500 SAMPLES/SEC„ 


FILTER 1 
CONSTANT GAINS 

RUN 2 

imput = noise 


KR 

RANGE 
VARIANCE 
REDUCTION RATIO 

KRR 

RANGE RATE 
VARIANCE 
REDUCTION RATIO 

T 

CONVERGENCE 
TIME (SEC.) 

tl = 0.7071 
K2 = 0.45 

6.276 X lO"^ 

2,605X10’^ 

23 

K1 =0.75 
K2 = 0.125 

7.192 xl0~^ 

2.023 X 10'® 

35 

K1 = 0.7 
K2 = 0.3769 

6.197 X lO'^ 

1.832 X 10'^ 

25 

K] = 0.49 
K2 = 0.98995 

2.493 X 10”^ 

1.250x10"® 

36 

K1 = G.75 
K2 = 0.45 

7.091 X 10‘^ 

2.614 x 10"^ 

23 

















TABLE V 


FILTER 1 RAW6E STANDARD DEV I AT IONS { FT. ) AT 2500 SAMPLES/SEC„ 


FILTER 1 
BLD 

RON 1 

i npuit = ramp 

RUN 2 

input = noise 

RUN 3 

input = ramp + noise 

INPUT 

RANGE 

STND. OEV. 

OUTPUT 
RANGE 
STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 
RANGE 
STND. DEV. 

' INPUT 
RANGE 

STND. DEV. 

OUTPUT 
RANGE 
STND. DEV. 

0.7778 
K2 = 0.3025 
b = 0.55 

0.1634 

0.1634 

110.7364 

0.9729 X 10*^ 

110.7339 

0.2251 

K1 = 0.7071 
K2 0.25 
b = 0.50 

0.1634 

0.1634 

110.7364 

.0.8841 X 10' V 

110.7339 

0.2468 

K1 = 0.6364 
K2 = 0.2025 
b = 0.45 

0.1634 

0.1634 


0.7951 X 10"^ 

110.7339 

0.2381 

K1 = 0.5657 
K2 = 0.160 
b = 0.40 

0.1634 

0.1632 

110.7364 

0.7056 X 10"^ 

110.7339 

0.2296 

K1 = 0.4950 
K2= 0.1225 
b = 0.35 

0.1634 

0.1633 

110.7364 

0.6157 X 10"^ 

110.7339 

0.2211 


r\> 

cn 

































TABLE VI 


FILTER 1 RANGE RATE STANDAPn DEVIATIONS(FT„/SEC. ) AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 

FILTER 1 

input = ramp 

input = noise 

input = ramp ■+* noise 


OUTPUT 

• OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 

BLD 

STND. DEV. 

STm. DEV. 

STNO. DEV. 

K1 = 0.7778 


■ 


K2 = 0.3025 

0.0 

0.3818 X 10"^ 

D.3818 X 10"^ 

b = 0.55 




K1 = 0.7071 




K2 - 0.25 

o,.o 

0.3155 X 10”^ 

0.3155 X 10"^ 





U = 0.6364 




K2 = 0.2025 

0.0 

0.2556 X 10‘^ 

0.2556 X 10'^ 

b = 0.45 




Kl = 0.5657 

i 



K2 = 0.160 

0.0 1 

0 -.2019 X 10"^ 

! 0.2019 X 10'^ 

b = 0.40 




Kl = 0.4050 




K2 = 0.1225 

0.0 

0.1544 X 10"^ 

0.1544 X 10"^ 

b = 0.35 



I 


ro 

















TABLE VII 


FILTER 1 MAXIMUM RANGE(FT.) AND RANGE RATE (FT. /SEC. ) ERRORS AT 2500 SAMPLES/SEG. 



RUN 1 j 

RUN 2 

RUN 3 


input = 

ramp 

input = 

noise 

input = ramp + noise 

FILTER 1 

MAX 

MAX 

MAX 

MAX 

MAX 

MAX 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

BLD 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

K1 = 0.7778 







K2 = 0.3025 

0.112 

0.901 X 1.0‘® 

■ 0.208 

0.246 

0.197 

0.246 

b = 0.55 







K1 = 0.7071 



* 

■ 

■ 



K2 = 0.25 

0.112 

0.128 X 10"® 

0.187 

0.223 

0.178 

0.223 

b = 0.50 



: 




K1 = 0.6364' 


■ ■ 





K2 = 0.2025 

0.112 X 10"^ 

0.161 X 10“^ 

0.172 

0.208 

0.161 

0.208 

b = 0.45 




■ 



K1 = 0.5657 







K2 = 0.160 

0.111 X 10“^ 

0.149 X lO-^ 

0.152 

0.199 

0.142 

0.199 

b = 0.40 







K1 = 0.4950 




1 



K2 = 0,1225 

0.119 X 10"^ 

0.757 X 10"'^ 

0.143 

0.188 i 

0.155 

0.188 

b = 0.35 
























FILTER 1 RANGE AND RANGE RATE 



VIII 


\NCE REDUCTION RATIOS AT 2500 SAMPLES/SEC. 


2 


= noise 


KRR 

RANGE RATE 
VARIANCE 
REOUCTTON RATIO 

T 

CONVERGENCE 
TIME (SEC.) 

1.189 X 10"^ 

24 

8.117 X 10‘® 

26 

5.328 X 10'^ 

29 

3.324 X 

32 

1.944 X 10“® 

J 

37 


ro 

VO 
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TABLE IX 


FILTER 1 RANGE STANDARD DEVIATIONS(FT, ) AT 6.25 SAMPLES/SEC. 



RUN 

1 

RUN 

2 

RUN 

3 


input = 

ramp . 

input = 

noi se 

Input = ramp + noise 1 

FILTER 1 

INPUT 

OUTPUT 

INPUT 

OUTPUT 

INPUT 

OUTPUT 


RANGE 

RANGE 

RANGE 

RANGE 

RANGE 

RANGE 

CONSTAflT GAINS 

STND. DEV. 

STND. DEV. 

STND. DEV. 

STND. DEV. 

STND. DEV. 

STND. DEV. 

K1 = 0.7071 

65.3068 

65.3067 

91.3215 

24.6606 

115.2178 

84.4420 

K2 = 0.45 





• 


K1 = 0.75 

65.3068 

65.3067 

91.3215 

19.8742 

115.2178 

78.3689 

K2 = 0.125 







K1 = 0.7 

65.30fi3 

65.3067 

91.3125 

23.6573 

115.2178 

83,4522 

K2 = 0.3769 

i 

! 






Kl = 0.49 

L 

65.3068 

65.3067 


24.2674 

115.2178 

82.9028 

K2 = 0.90995 




i 



ICl = 0.75 

65.3068 

65.3067 

91.3215 

24.8739 

115.2178 

84.1560 

K2 = 0.45 

■ 1 








ro 






























TABLE X 


FILTER 1 RANGE RATE STANDARD DEVTATIONS(FT./SEC. ) AT 6.2o SAMPLES/SEe. 


riLTl£R 1 
CONSTANT GAINS 

RON 1 

input = range 
OUifPUT 
RANGE RATE 
STND, DEV. 

RUN 2 

iinpuit = noise 
OUTPUT 
RANGE RATE 
STND. DEV. 

RUN 3 

i!nput = ramp * noise 
OUTPUT 
RANGE RATE 
STND, DEV. 

KI = G.7071 




K2 = 0.4'5 

0.0 

11.4052 

11.4052 

Kl = 0.75 




K2 = 0.125 

0.0 

3,2672 

3.2671 

Kl = 0.7 

0,0 

9.8217 

9.8218 

K2 = 0.3769 




K 1 = 0,49 

'g.o 

24.8694 

24.8694 

K2= 0.98995 




Kl=0.75 




K2=0.45 

O.G 

11.2231 

11,2231 


CO 

















TABLE XI 


FILTER 1 MAXIMUM RANGI(FT.) AND RANGE RATE (FT. /SIC. ) ERRORS AT 6.25 SAMPLES /SEC. 



RON 1 

RUN 2 

RUN 3 


Input - 

ranm 

input = 

noise 

Input = ramp + noise 

F IiLTER 1 

MAX 

MAX 

MAX 

MAX 

MAX 

MAX 

RANGE 

RA'NGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

G0N5TANT GAIKS 

ERROR 

ERROR 

ERROR 

ERR0R 

ERROR 

ERROR 

K1 = 0.7G71 

4 .4000 

0.102 X 10'^ 

62.8718 

28.0303 

57.3518 

28.0303 

K2 = 0.45 

■ 

; 

■ 





Kl=-0.75 

4.4800 

0.126 X 10"^^ 

44.8687 

5.7304 

49.3487 

5. 7304 

K2 = 0.125 







K1 =0.7 

4.4800 

0.116 X 10'^^ 

51.5882 

22.9018 

56.0682 

22.9018 

K2 = 0.3769 







K1 = 0.49 

•- 

0.483 X H0~^^ 

1 



*■ 



4.4800 

5.7.2146 

63.0950 

61.6946 

63.0950 

K2 = 0.98995 







K1 = 0.75 

4.4800 

G.112 X 10"^^ 

53.3341 . i 

27.4471 

57.8141 

i 

27.4471 

X2 = 0.45 








-fi 

4 ^ 

































TABLE XII 


FILTER 1 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SECOND 


FILTER 

COMSTAHT GAINS 

RUN 2 j 

input = noise 


KR 

RANGE 
VARIANCE 
REDUCTION RATIO 

KRR 

RANGE RATE 
VARIANCE 
REDUCTION RATIO 

T 

CONVERGENCE 
TIME (SEC.) 

K1 = 0.7071 
K2 = 0.<1'5 

7.292 X 10'^ 

1,560 X lO’^ 

23 

K1 = 0.75 
K2 = 0.125 

4.736 X 10"^ 

1.280 X 10'^ 

36 

Kl = 0-7 
K2 = G.3769 

6.711 X lO"^ 

1.157 X 10‘^ 

25 

Kl = 0.4'9 
K2 = 0.90995 

7.062 X 10‘^ 

7.416 X 10'^ 

38 

Kl = 0.75 
K2 ^ 0.45 

7.41;9 X 10"^ 

1.510 X 10"^ 

23 










TABIE XIII 


FILTER 1 range STANDARD DEVIATIONS ( FT„ ) AT 6.25 SAMPLES/SEC. 


FILTER 1 
BED 

RUM 1 

input = ranip 

RUN 2 

input = no1S6 

RUN 3 

i'HPut = ramp + noise 

INPUT 
RANGE 
STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND, DEV. 

OUTPUT 
RANGE 
STND. DEV. 

INPUT 
RANGE 
STND. OEV. 

OUTPUT 

RANGE 

STND. DEV. 

K1 = 0.7778 
K2 = 0.3025 
b = G.55 

65.3068 

65.3067 

91.3215 

23.1994 

115.2178 

82.0274 

K1 = 0.7071 
K2 = 0,25 
b = 0.50 

65.3068 

65,. 3067 

■ 

91.3215 

21.6291 

115.2178 

81.0868 

Kl = 0.6354 
K2 = 0.2025 

65.3068 

65.3067 

91.3215 

19.9024 

115.2178 

79.8818 

K1 - 0.5657 
K2 = 0.160 

b =0.40 ^ 

65.3068 

65.3067 

91.3215 

17.9741 

115.2178 

78.3595 

Kl = 0.4950 1 
K2 = 0.1225 1 
b = 0.35 

^ . ■■ ■ J 

65.3068 

i 

65.3067 

91 ! 32 15 

15.7742 

115.2178 

76.4249 


cr> 
























TABLE XIV 


FILTER 1 RANGE RATE STANDARD DEVIATIO:NS{FT./SEC, ) AT 6.2'5 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


input = ramp 

Input = noise 

input = ramp + noise 

FILTER 1 

OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE PATE 

BLf) 

STND.. 0EV. 

STNB. DEV. 

STND. DEV. 

K.1 = 0.7778 




K2 = 0.3025 

0.0 

7.7822 

7.7822 

b = 0.55 




K1 = 0.7071 




:K2 = 0.25 

0.0 

6.7482 

6.7482 

-b = 0.50: 




K1 = 0.6364 




K2 = 0.2025 

0.0 

5.7662 

5.7662 

b = 0.45 




Ki = G.5657 




K2 = 0.160 

0.0 

4.8206 

4.8206 

b = 0.40 




Kl = 0.4950 




K2 = 0.1225 

0.0 

3.8874 

3.8874 

b = 0.35 
















TABLE XV 


FILTER 1 MAXIMUM RANGE ( FT. ) AND RANGE RATE(FT./SEC.) ERRORS AT 6.25 SAMPLES/ SEC. 



RUIN 

1 

RON 2 1 

RUIN 

3 


iinpuit = 

^'anip 

input - 

noise 

input = ranip + noise 

FILTER 1 

MAX 

rax 

MAX 

MAX 

■I 

MAX 

MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

BLB 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

XI = 0.7778 







K2 = 0.3025 

4.4'8D 

0.209 X 10"^^ : 

• 50.988 

16.971 

55,468 

16.971 

b =0.55 

. 






K1 = 0.7071 







K2 = 0.05 

4.480 

0..245 X 10'^^ 

48.114 

13.780 

52.594 

13.780 

b = 0.50 







XI = 0.6364 







K2 = 0.2125 

4,480 

0.270 X 

44.710 

10.906 

49.190 

10.906 

b = 0.45 







XI = 0.5657 

1 







1 

K2 = 0.160 
b = 0.40 

4.48G 

0.329 X lO"^^ 

40.638 

8.570 

45.168 

8.570 

XI = 0.4'950 

' ' 






K2 =0.1225 

4,480 

Q.330 X 10"^^ 

36.010 

6.730 

40.490 

6.730 

b =0.35 








-P* 

00 





















TABLE XVI 


FILTER 1 RANGE AND RAWGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEC. 


FILTER 1 

BID 

RUN 2 

iinpuit = noise 


KR 

RANGE 
VARIANCE 
REIDCTION RATIO 

KRR 

RANGE RATE 
VARIANCE 
REDbGTION RATIO 

T 

C0NVERGENCE 
TIME (SEC.) 

XI = 0.777a 
K2 = 0.3Q2 d 
b = 0.55 

6.454 X 10"^ 

7.262 X 10'^ 

23 

XI = 0.7071 
K2 = 0.25 
b = 0.50 

5.610 X 10'^ 

5.460 X 10"^ 

26 

K1 = 0.5364 
K2 = 0.2025 
5 =0.45 

4.750 X 10"^ 

3.987 X 10"^ 

28 

XI = 0.5657 
K2 = 0.160 
b = 0.40 

3.874 X 10"^ 

2.770 X 10"^ 

31 • 

K1 = 0.4950 
K2 = 0.1225 
b = 0.35 

2.984 X liO"^ 

1.812 X 10"^ 

36 
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RUN 1. With the ramp input, the recorded range error term repre- 
sents the change in range over one sampling period. This can be 
corrected by multiplying the numerator of the transfer function 
Hj^ (equation (3)) by the delay term This has the effect 

of shifting the system back in time one sampling period. The 
resulting transfer function is 

Li „ RE - z“MK1DT + K20T^) - z“2(K1DT) 

"r RM 1 + z-^X-1 "+ KIDT + K2D;^) + - KIDT}’ 


2, Filter 2 

Simulation data for Filter 2 is presented in Tables 
XVII - XX for 2500 samples/second and in Tables XXI - XXIV for 
the 6.25 sampTes/second sampling rate, A graphica') display of 
the data is shown for each set of parameter pairs a and 3 at 
2500 samples/second in Figures 31-44, A graphical display of the 
data for the 6.25 samples/second sampling rate is shown in Figures 
45-53, Graphs for two sets of parameter pairs are omitted because 
they yield systems whose convergence times exceed the duration 
of the simulation run. 

The alpha-beta filter is simulated using six sets of 
a ,6 pairs that Tie on the curve described by (10) within the 
region of stability defined by: 


a > 0, g > 0 and 2a + g<4. 


( 22 ) 



Near the upper end of the curve;at the point a = 1.0 and & - 1.0, . 
the filter is characterized by fast maneuver-following ability 
but poor noise rejection. Conversely, near the Tower end of the 
curve, as a and S approach zero, the filter is characterized by 
improved noise rejection at the expense of maneuver- foil owing 
ability. 

For the 2500 sampWsecond tests, relative manuever-fol low- 
ing capabilities can be determined by examining the convergence 
time. This is the time required by the filter to recover from 
the transient at the start of the simulation. For larger values 
of {such as a = 1.0), convergence time is negligible. However, 
the range variance reduction ratio (KR) is one, meaning there is 
no noise attenuation. As a and g are decreased, convergence time 
increases indicating the inertia of the system is increasing and 
more time is required to recover from transients. While this 
presents a reduction in maneuver-fon owing capabilities it also 
means the system is becoming more insehsitive to the higher fre- 
quency components of noise. Hence, there is a corresponding 
decrease in range and range rate variance reduction ratios. 

These observations are even more readily seen in the 
6.25 sample^^seeond simulation runs. The reduction in sampling 
rate causes an increase in convergenee time. As an example, for 
a ?= 0.05 and g = 1.282 x 10”^, the convergence time at 2500 
samples/second is less than one second. However, at 6.25 samples/ 
second that figure is approximately 56 seconds. In fact, the last 
three entries in the tables for the 6.25 samples/second tests all 
represent a, 6 pairs whose convergenee times exceed the 192 



TABLE XVII 


FILTER 2 RANGE STANDARD DEVIATIONS ( FT. ) AT 2500 SAMPLES/ SEC. 


FILTER 2 
alpha - beta 

RUN .1 

Inout = Fenm_ 

RUN 2 

input = noise 

RUN 

To out = ram 

z 

p + noise 

INPUT 

RANGE 

STND. BEV. 

r 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

QUTPIOT 

RANGE 

STND. DEV. 

INPUT 
RANGE 
STND. DEV, 

OUTPUT 

RANGE 

STND, DEV. 

a = 1.0 

a = 1.0 

0.1634 

0.1634 

110.7364 

• 

110.7364 

110.7339 

110.7339 

a = 0.2 
3 - 0.022222 

0.1634 

0.1634 

110.7364 

29.0844 

110.7339 

29.0377 

a = 0.G5 

3 = 0.1282 X lO"^ 

0,1634 

0.1634 

110.7364 

7.5559 

110.7339 

7-6332 

■ 1 

a = O.Q01 
3 = 5.0025 X 10'^ 

0.1634 

1 

0.1634 

110.7364 

i 

0.3205 

110.7339 

0.4751 

ft = 0.00030 
’ 3 = 4.50067 X 10‘® 

O’. 1634 

0.1634 

11Q..7364 

1 

0.9379 X 

110.7339 

0.2518 

a = 0.OQ021 
3 = 2.2G523 x 10”® 

0.1634 

0.1632 

110.7364 

0.6538 X 10^ 

-110.7339 

0.225 




TABLE XVIII 

FILTER 2 RANGE RATE STANDARD DEVrATIONS(FT./SEC . ) AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


input = rainp 

Input = noise 

input = ramp + noise 

FILTER 2 

OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 

alpha - beta 

STND. 0EV. 

STND. DEV. 

STND. DEV. 

a = l.G 
B = 1.0 

0.0 

412275.80 

412275.81 

11 

p 

- 



0 = 6.0222222 

0.0 

8741.49 

8741.49 

a ~ 0*05 

■ 



B = G.1282 X 10"^ 

0.0 

4.45.064 

445.0637 

a = G.001 
0 = 5.0025 X 10“^ 

0.0 

0.3965 

0.3965 

a = 0.00030 


0.3550 X 10“^ 

0.355 X 10^^ 

0 = 4.50057 X 10"® 

0.0 

a = 0,00021 




0 = 2.20523 X 10"® 

0.0 ' 

0.1739 X 10"^ 

0.1739 X 10"^ 


er* 

tn 



















TABLE XIX 


FILTER 2 MAXIMUM RA'NGE{FT.) AND RANGE RATE (FT. /SEC. ) ERRORS AT 2500 SAMPLES/SEG. 


FILTER 2 
alpha - beta 

RON 1 

input = ramp 

RUN 2 

input = noise 

RUN 3 

input = ramp + noise 


mx 

RANGE RATE 
ERROR 

MAX 

RANGE 

ERROR 

MAX 

RANGE RATE 
ERROR 

MAX 
RANGE 
ERROR ■ 

MAX 

RANGE RATE 
ERROR 

a = 1.0 

e = 1.0 

o.o’ 

0.578 X 10'^ 

. 

254.932 

989103.5 

254.932 

989103.6 

« = 0.2 
S = 0.022222 

0.0 

0.226 X 10"^® 

96.237 

21226.7 

96.237 

21226.7 

ct = 0.05 

0=0.1282x10"^ 

0.909 X 10"^^ 

0.476 X 10"® 

37.824 

1792.05 

37.'824 

1792.05 

■ ' 1 

a = 0.001 1 

3 = 5.0025 X 1Q“^ 

0.546 X 10" 

’ . ■ .! 

0.166 X 10“® 

0.780 

, 1 

m 

0:780 

1.409 

a = G. 00030 
3 = 4.50067 X 10"® 

0 .425 X 10"^ 

0.414 X 10"® 

0.199 

0.236 

0.199 

0.236 

a = 0.00021 
3 = 2.20523 X 10"® 

0.134 X 10"'^ 

0.670 X 10"'* 

0.133 

0.194 

0.133 

i 

0.194 

































TABLE XX 


FILTER 2 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 2500 SAMPLES/SEG. 



RUN 

2 

PREDICTED BY 



input = 

noise 

BENE.0IGT AND BORDNER 


FILTER 2 
alpha - beta 

KR 

RANGE VARlANeE 
REDUCTION RATIO 

KRR 

TANGE RATE VARIANCE 
REDUCTION RATIO 

KR 

KRR 

T 

CONVERGENCE 
TIME (SEC.) 

a = 1.0 
& = 1.0 

l.Q 

1.386 X 10^ 

1.0 

1.25 X 10^ 

<1 

a = 0.2 
3 = G. 022222 

6.898 X 10"^ 

' . - - . 

6.231 X 10^ 

1.534 X 10"^ 

8.521 X 10^ 

<1 

a = G.05 

e = 0.1282 X, 10'^ 

4. 656 X 10"^ 

1.615 X 10^ 

3.782 X 10"^ 

1.056 X 10^ 

<1 

a = 0.001 
8 = 5.0025 X 10'^ 

8.377 X 10"® 

1.282 X I'O"® 

7.501 X 10"^ 

7.824 X 10"^ 

8 

a = 0.00030 
8 = 4.50067 X 10‘® 

7.174 X 10"^ 

1..028 X 10"^ 

2.250 X 10"^ 

2.110 X 10"® 

24 

a = 0.00021 
8 = 2.20523 X 10^® 

3.486 X 10'^ 

2.466. X 10"® 

1.575 X 10"^ 

7.237 X 10"® 

35 
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FIGURE 36 - Filter 2 at 2500 samples/see. 
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TABLE XXI 


FILTER 2 RANGE STANDARD DEVlATr0NS(FT. ) AT 6.25 SAMPLES /SEC. 



RON 1 

RUN 2 


RUN 3 


_ _ _ Input = 

^ raniD 

iinout = 

noise 

inout = 

ramp + noise 

FILTEK 2 
alpha - beta 

INPUT 
RANGE 
STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 
RANGE 
STND. DEV 


OUTPUT 

RANGE 

STND. DEV. 

a = l.O 

e = 1.0 

65.3068 

65.3068 

91.3215 

91.3215 

115.2178 

115 .2178 

ee = 0.2 
$ = 0.0222.22 

65.3068 

65.3068 

91.3215 

30.0503 

115.2178 

83.4223 

a = 0*05 

|3 = 0.1282 X 10'^ 

65.3058 

65.3068 

i 

91.3215 

8.7326 

115.2178 

69.3696 

a = 0.001 








0 = 5.0025 X lO"^ 

65.3068 

i 

91.3215 

0.2055 

115.2178 



— — — — — rj 

a = 0.00030 

65.3068 

i 

91.3215 

0.6791x10'^ 

115.2178 

' 

0 = 4.50067 X 10"® 

i 

j 

a = 0.00021 

65 .3068 


91.3215 

0.4817x10'^ 

■ 

115.2178 


3 = 2.20523 X 10‘® 
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TABLE XXIII 


\ ■ 
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FILTER 2 MAXIMUM RANGE{FTj AN0 RANGE RATE {FT./SEC.) ERRORS AT 6*25 SAMPLES /SEC* 



RON 

1 

RUN 2 


RUN 3 


input = 

ramp 

iinput = 

noise 

input “ 

ramp + noise 

FILTER 2 

MAX 

MAX 

MAX 

MAX 

MAX 


MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 


RANGE BATE 

alpha - beta 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 


ERROR 

a = 1.0 

o!o 

0.632x10"^^ 

230.565 

1613.121 

230.565 


1613.121 

P = 1.0 








a = 0.2 

0.0 

0.233x10"^^ 

71.466 

51,347 

71.466 


51.347 

0 = 0.Q22222 








a = G,05 

p 

0.279x10"® 






0 = 0.1282 X 10'^ 

0. 646x10"“^ 

. 23.015 

3.665 

23.015 


3.665 

a = 0.001 



5.503 

0.136x10"^ 




e = 5.0025 X 10"^ 

i 





a = O*OQ030 ' 

‘ 


j 

1.964 





e = 4.50067 X 10"® 



0.172x10 





a = 0.00021 



; 

0.879x10"® 




Q 



1.404 




3 = 2.20523 X 10"^ 









GO 
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7:able XXIV 


FILTER 2 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEC. 



RUN 2 1 

PREDICTED BY 



input = 

noise 

BENEDICT AND BORDNER 


FILTER 2 
alpha - beta 

KR 

RANGE VARIANCE 
REDUCTION RATIO 

KRR 

lANGE RATE VARIANCE 
REDUCTION RATIO 

KR 

KRR 

T 

CONVERGENCE 
TIME (SEC.) 

a = 1.0 
3 = 1.0 

1.0 

67.778 

1.0 

78.125 

<1 

a = 0.2 
3 = 0.022222 

1.0B3 X 10"^ 

4.345 X lO'^ 

1.554 X 1§‘^ 

5.392 X 10'^ 

15 

ri ~ 0.05 

3 = 0.1282 X 10‘^ 

9.144 X 10'® 

2.842 X 10'^ 

3.782 X 10'^ 

6.587 X 10"'’ 

56 

a = 0.001 

•3 -- 5.002B X 10"^ 

5.064 X 10‘® 

4.160 X 10"^^ 

7.501 X lO"'* 

4.890 X 10"® 

>192 

a = 0.00030 
3 = 4.50067 X 10‘® 

5.530 X 10"^ 

4.497 X 10“^® 

1 

1 

-4 

2.250 X 10 ^ i 

1.319 X 10"’® 

>192 

a = 0.00021 
3 = 2.20523 X 10'® 

2.782 X 10~'^ 

1.132 X 10“^® 

1.575 X 10"^ 

4,523 X 10"^^ 

>192 
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second duration of the tests. For this reason steady state values 
for output range and range rate standard deviations are not given. 

The range variance reduction ratio is not a function of 
sampling rate as can be seen by comparing that value for both 
rates. The range rate variance reduction ratio, however, is 
proportional to the square of the sampling rate and demonstrates 
a decrease with a reduction of sampling rate. 

Overall, data from the simulation runs for the alpha- 
beta filter agrees with the variance reduction ratios predicted 

5 

by Benedict and Bordner. Deviations become apparent as alpha 
and especially beta decrease in magnitude^ 

3. Filter 3 

Simulation data for Filter 3 is presented in Tables 
XXV - XXVIII for 2500 samples/second and Tables XXIX - XXXII 
for the 6.25 samples/second sampling rate» A graphical display 
of this data is shown in Figures 54 and 55 for 2500 samples/second 
and Figures 56 and 57 for 6.25 samples/second. 

The N sample smoother employed to derive range in Filter 
3 smooths over sixteen samples. That is a sufficient number of 
samples to provide good filtering, However, its performance is 
not readily apparent from the data. Since the smoother provides 
a range value that lags the present sample instant by 7,5 sampling 
periods, updating is necessary. To accomplish this a correction 
term is determined from the estimated range rate calculated for 
that sample instant. This updating introduces error into the esti- 
mated range. 
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For simulation RUN 2, since both the range and range rate 
segments of Filter 3 are linear systems and the input is a random 
signal, the estimated range before correction for the 7,5 sample 
period lag (RE lagging) and the estimated range rate (ROOT) are 
random signals. The operation of correction for lag can then be 
described in terms of variances as 


^RE lagging * ^RDOT^^ ^RE 


( 22 ) 


where K is the constant 7„5 * DT. Using the values for estimated 


range standard deviation (<?[^g) and estimated range rate standard 


deviation {<7|^[)0j) listed under RUN 2 in the tables for Filter 3, 
the standard deviation of the lagging estimated range (aj^^ 
can be obtained. This value is 3. 9291m. (12.8908 ft.) and 
7.0905m, (23.2627 ft.) for the 2500 and 6.25 samples/second simula- 
tions, respectively. The associated range variance reduction ratios 


are 


’^'^2500 


_/ 12.8908s 2 
■^110.7364^ 


= 1.355 X 10'^ 


(23) 


1^6.25 = 'Irilf'' = 

This is in agreement with the variance reduction ratio Of an 
N sample smoother. A eomparisOn can be made of these ratios 
to those given in the tables to demonstrate the effect the up- 
dating procedure has on noise rejection in range detemti nation. 



! 
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TABLE XXV 


FILTER 3 RANGE STANDARD DEVIATIONS(FT. ) AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


input = 

ramp 

input = 

noise 

input = ramp + noise 


INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 
RANGE 
STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

FILTER 3 

0,1634 

0.1634 

110.7364 

34.7179 

t — : 

110.7339 

34.6758 


TABLE XXVI 


Fin 

■ 

ER 3 RANGE RATE STANDARD C 

RUN 1 
•OUTPUT 
RANGE RATE 
STND. DEV. 

iEVIATIONS(FT./SEC.) AT 251 

RUN 2 
OUTPUT 
RANGE PATE 
STND. DEV. 

DO SAMPLES/SEG. 

RUN 3 
OUTPUT 
RANGE RATE 
STND. DEV. 

FILTER 3 

0.0 

10745.3359 

10745.3359 


TABLE XXVII 


FILTER 3 MAXIMUM RANGEfFT.) AND RANGL RATE(FT./SEG.) ERRORS AT 2500 SAMPLES/SEC. 



RLIM 

1 

i. 

RUN 2 

c: 


input = 

ramp 

input = 

noise 

input = ramp + noise | 


MAX 

MAX 

MAX 

MAX 

I4AX 

rnAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 


ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

FILTER 3 

0.364 X 10"^^ 

0;207 X 10"® 

114.486 

28399.3 

114.4859 

28399.3 


TABLE XXVII! 

FILTER 3 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 2500 SAMPLES/SEC. 



RUN 2 



KR 

KRR 

T 


range variance 

range rate variance 

Convergence 


reduction ratio 

reduction ratio 

Time (Sec.) 

FILTER 3 

9.829 X 10"^ 

9.416 X 10^ 

<1 



169. QD 














TABLE XXIX 


FILTER 3 MNGE STANDARD DEVIATIONS(FT. ) AT 6.25 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


i’nput = 

ramp 

input = 

noi se 

imput = rainp + noise 


INPUT 

RANGE 

SEND. DEV. 

OUTPUT 

RANGE 

SEND. DEV. 

■ 

INPUT 

RANGE 

SEND. DEV. 

OUTPUT 
RANGE 
STND. DEV. 

INPUT 

RANGE 

SliND. DEV. 

OUTPUT 

RANGE 

STNO. DEV. 

FILTER 3 

65.3066 

65.3068 

91.3215 

36.3516 

115.2178 

86.9980 


TABLE XXX 


FILTER 3 RANGE RATE STANDARD DEVIATIONS(FT,/SEG, ) AT 6.25 SAMPLES/SEG. 



RUN 1 

RUN 2 

RUN 3 


•OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 


STND. DEV. 

STND. DEV. 

STND. DEV. 

FILTER 3 

0.0 

1 

23.2780 

23.2700 

i 




TABLE XXXI 

FILTER 3 MAXIMUM RANGE (FT. ) AND RANGE RATE(FT./SEC. ) ERRORS AT 6.25 SAMPLES/SEG. 



TABLE XXXII 


FliTER 3 


RANGE RATE VARIANGE REDNCTION RATIOS AT 6.25 SAMPLES/ SEC. 
RUN 2 


KR 

range variance 
reduction raitio 


KRR 

range rate variance 
reduction ratio 


Convergence 
Time (Sec.) 


FILTER 3 


1.584 X 10 


-1 


6.497 X 10 


rZ 


. 3 
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The range rate detertn’ini'ng section of Filter 3 Is a'cas- 
caded simple average smoother. The range rate variance reduction 

g 

ratio is given by 


RATE OUTPUT _ ^ h ^ 

°^N6E INPUT Tnpig^TP' * 

where is the number of samples over which the first differences 
are taken (n^ = 12) and r \2 is the number of first differences 
that are smoothed (ng = 6)» Simulation runs at both 2500 and 6.25 
samples/second support the above relation. Operation at 6.25 
samples/seGond gives improved results due to the dependence of 
the variance reduction ratio on the square of- the sampling period 
DT. 


4. Filter 4 

Simulation data for Filter 4 is presented in Tables 
XXXI I I - XXXVl for 2500 samples/second and in Tables XXXVII - 
XL for 6,25 samples/second. A graphical display of this data 
is shown in Figures 58 and 59 for 2500 samples/second and in 
Figure 60 and 61 for 6.25 samples/second. 

The digital tracking filter smoothing over four Samples 
has the transfer functions 


_ RE _ 1 + + z"^ + z"^ 

R " RM ■ 1 % 0,75 + O.ZSz-^ + tr.25z-^ + 0.251"^ (26) 


and 


TABLE XXXIII 


FILTER 4 RANGE STANDARD DEVIATIONS ( FT, ) AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 

2 

RUN 3 


Input = 

ramp 

Input = 

noise 

Input = ramp + noise 


INPUT 

OUTPUT 

INPUT 

OUTPUT 

INPUT 

OUTPUT 


RANGE 

RANGE 

RANGE 

RANGE 

RANGE 

RANGE 


STKO. BEV. ; 

STNO. DEV. 

STND. BEV. 

STNO. DEV. 

STND. DEV. 

STND. DEV. 

FILTER 4 

0.1634 

0.1634 

110.7364 

144. 8444 

110.7339 

144.8549 


TABLE XXXIV 


FILTER 4 RANGE RATE STANDARD DEVIATIONS{FT./SEGo) AT 2500 SAMPLES/SEG. 



RUN 1 

RUN 2 

RUN 3 


OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 


STND. DEV. 

STNO. DEV. 

STNO. DEV. 

FILTER 4 

0.0 

247907.688 

247907.683 















J -Jf 

■'■i 


1 


I 


ft: ■■ 


t 




TABLE XXXV 

fllTER 4 MAXIMUM RANGE { FT „) AND RANGE RATE(FT./S£G. ) ERRORS AT 2500 SAMPLES/SEG. 



RUN 

1 

RUM 2 

RUN 3 


input = 

ramp 

input = 

noise 

input “ ramp + noise 


MAX 

MAX 

MAX 

MAX 

MAX 

MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 


ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

FILTER 4 

0.182 X 10'^^ 

0,320 X 10"^ 

290.7632 

501279.4 

290.7632 

501279.375 


TABLE XXXVI 

FILTER 4 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 2500 SAMPLES/SEC. 



RON 2 



KR 

KRR 

T 


range variance , 

range rate variance 

Convergence 


reduction ratio 

reduction ratio 

Time (Sec.) 

FILTER 4 

1.712 

5.012 X 10® 

<1 








10«00 


FIGURE 59 - Fil 







TABLE XXXVII 


FILTER 4 RANGE STANDARD DEVIATIONS(FT. ) AT 6.25 SAMPLES/SEC. 



RUN 1 

HUN 2 

RUN 3 


iinput = 

ramp 

iinput = 

noise 

input = ramp + noise 


||H|| 

OUTPUT 

RANGE 

STND. DEV. 


OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 


FILTER 4 

65.3068 

65.3068 

91.3215 

15? ,1930 

1 

L_ 

115.2178 

178.7200 


TABLE XXXVIII 


FILTER 4 RANGE RATE STANDARD DEVIATIQNS(FT./SECo ) AT 6.25 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 


STND. DEV. 

STND. DEV. 

STND. DEV. 

FILTER 4 

0.0 

693.4480 

693.4480 















TABLE XXXIX 

FILTER 4 MAXIMUM RANGE{FT.) AND RANGE RATE(FT./SEG. ) ERRORS AT 6.25 SAMPLES/SEC. 



RUN 

1 

RUN 2 

RUN 3 


iinput - 

ramp 

Input = 

noise 

input = ramp + noise 


MAX 

MAX 

MAX 

MAX 

MAX 

MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RMiGE 

RWGE RATE 


ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

FILTER 4 

0.182 X 10"^^ 

0.580 X 10’^^ 

323.3982 

1553.83 

323.3982 

1553.82 


TABLE XL ' 

FILTER 4 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEC. 



RUN 2 



KR 

KRR 

T 


range variance 

range rate variance 

Convergence 


redut ^i:on ratio 

reduction ratio 

Time (Sec.) 

FILTER 4 

2.9603 

57.6609 

[ 

2 
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,, _RDOT _ 1-0.75z‘*-1.5z”2-2.25z”5-4z"'*-2.25z’5-1.5z“®-0.75z‘”^ /o-rv 
^RR''W " 4DTll-0.75z-^ + 0.25z-^+ 0.25z-’ + 0.25r^) ' 


The variance reduction ratios, KR and KRR, can be determined from 
the weighting sequences h|^(i) and hpp(i) by 


K= h"(1), 


(28) 


The program presented in Appendix D is used to evaluate this 
series* The resulting variance reduction ratios calculated are 

KR = 48 (29) 

KRR(250O samples/sec.) = 1.47 x 10® (30) 

KRR(6.25 samples/sec. ) = 916. (31) 

The variance reduction ratios presented in the tables for Filter 4 
are ealculated from actual simulation data. Although these are 
slightly better results than determined above, smoothing over 
four samples does not offer an adequate level of noise rejection. 

Filter 5 

Simulation data for Filter 5 is presented in Tables 
XLI “ XL IV for 2500 samples/seeond and in Tables XLV - XLVIII 
for 6.25 samples/ Second. A graphical display of this data is 
presented in Figures 62 and 63 and in Figures 64 and 6:5 for 6.25 
samples/seeond. 


115 


Filter 5 is a digital tracking filter employing skip 
sampling over sixteen samples. It is characterized by the transfer 
functions 

. RE _ 23.5 + 23.5z‘^+...+23.5z"^5-7.5z"'® -...-7.5z'^^ 

’ RH 256 (32) 


and 


„ „ ROOT _ 1 + (33) 

”RR ■ lifr " “ 256 DT 

The program to evaluate variance reduction ratios (listed in 
Appendix D) using 




detemined the ratios: 


KR = 1.72 X 10"^ (34) 

KRR(2§00 samples/seG. ) 5,74 x 10^ (35) 

KRR(6,25 samples/sec,) = 3.59 x 10“^ (36) 

These values are in close agreement with those derived from 
simulation data and presented in the tables for Filter 5. As 
with this type of filter the range rate variance reduction ratio 
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Is proportional to the square of the sampling rate while the range 
variance reduction ratio is indep^ident of sampling rate. 
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TA'BLE XL I 


FILTER 5 RANGE STANDARD DEVIATI0NS(FT. ) AT 2S00 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 3 


Input = 

= ramp 

Input = 

noise 

input = ramp + noise 


INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. . 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

: STND. DEV. 

FliLTER 5 

2.6123 

2.6123 

106.1125 

39.3340 

1G6.1860 

39.3009 


TABLE Kill 


FliLTER 5 RANGE RATE STANDARD DEVIATI0NS(FT./SEG.) AT 2500 SAMPLES/SEG. 


• 

RUN 1 

RUN 2 

RUN 3 


OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 


STND. DEV. 

STND. DEV. 

STND. DEV, 

FILTER 5 

0.0 

5693.6289 

5693.6289 
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TABLE XLIII 

FILTER 5 MAXIMUM RANGE(FT.) AND RANGE RATE(FT./SEC. ) ERRORS AT 2500 SAMPLES/SEC. 



RUN 1 

RUN 2 

RUN 

3 


input = 

ramp 

input = 

noise 

Input = ,rainp + noise I 


MAX 

MAX 

MAX 

MAX 

MAX 

MAX 


RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 

RANGE Smi 


ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

FILTER 5 

0 0.0 

O;.207=x 10"® 

122.399 

1 

15702.340 

122.399 

15702.340 


TABLE XL IV 

FILTER 5 RANGE AND RANGE RATE VA'RIANGE REDUCTION RATIOS AT 2500 SAMPLES/SEe. 



RUN 2 



K'R 

KRR 

T 


range variance 

range rate variance 

Convergence 


reduction ratio 

reduction ration 

Time (Sec.) 

FILTER 5 

1.374 X 10‘^ 

2.879 X 10^ 

<1 
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FIGURE 63 - Filter 5 at 2500 samples/seG. 
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TABLE XLV 

FILTER 5 RA'NGE STANDARD DEVIATIONS (FT. ) AT 6.25 SAMPLES/SEC. 



RUN J 

RUN 2 

RUN 3 


input = 

' ramp 

input - 

noise 

input = ramp + noise 


INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STNB. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. 

flOER 5 

1D44.906 

1044.96 

112.3026 

41.1501 

1083.432 

1044.917 


TABLE XLVI 

FILTER 5 RANGE RATE STANDARD DEVIATI0NS(FT./SEG. ) AT 6,25 SAMPLES/SEG. 



RUN 1 

RUN 2 

RUN 3 


OUTPUT 

OUTPUT 

OUTPUT 


RANGE RATE 

RANGE RATE 

RANGE RATE 


STND. DEV. 

STND. DEV. 

STND. DEV. 

FILTER 5 

0.0 

15.0085 

15.0085 



TABLE XiVII 


FILTER 5 MAXIMUM RANGE (FT.) AND RANGE RATE (FT. /SEC. ) ERRORS AT 6.25 SAMPLES/SEC. 



RUf 

1 

RUN 2 


RUN 3 


input = 

ramp 

input = 

noise 

input = 

.pamp + noise 


MAX 

MAX 

MAX 

MAX 

MAX 


MAX 

■ 

RANGE 

RANGE RATE 

RANGE 

RANGE RATE 

RANGE 


RANGE RATE 


ERROR 

ERROR 

ERROR 

ERROR 

ERROR 


ERROR 

FILTER 5 

0.9G9 X 10"^^ 

0.278 X 10"^^ 

122.002 

37.450 

122.0018 

37.450 


TABLE XLVIII 

FILTER 5 RANGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEC. 



RUN 2 



KR 

range variance 
reduction ratio 

KRR 

range rate variance 
reduction ration 

T 

Convergence 
Time (Sec.) 

FILTER 5 

1.342 X 10"^ 

1.786 X W~^ 

5 


ro 

ro 


FIGIilRI 64 - Filter 5 at 6.25 samples /sec. 












IV. MICROPROCESSOR SIMULATIONS 


Two of the filter designs are selected to be implemented on a 
microprocessor based system. They are Filter 3 and Filter 5. Simula- 
tions are conducted in the same manner as is done on the IBM 370. The 
same range input signals are employed, and statistical analysis of the 
outputs is performed as was done in the IBM 370 simulation runs. This 
allows direct comparisons to be drawn. 

the simulations are not conducted in real time. Considerations 
are taken for execution time, however, and are presented in the test 
results. 


A. Description of Micropre c essor Systein 

The mfcroprocessor simulations are conducted on an expanded 
version of the Motorola MEK 68O0D2 Evaluation Kit 11* a microcomputer 
system based on the Motorola MC 6800 microprocessing unit. The details 
of the revisions made to the system are outlined in the Motorola appli- 

7 

cation note "MEK 6800D2 Microcomputer Kit System Expansion Techniques". 
The basic system memory of 512 bytes of static RAM is replaced with 32 
kilobytes of dynamic RAM. The MINIBUG III monitor ROM and its associated 
serial interface are added co-resident to the basic system's JBUG monitor 
ROM and Serial interface. Under MINIBUG control the system is connected 
via RS^232 interface to a data terminal . Under JBUG control the audio 
cassette interface provides an inexpensive means of permanent storage. 
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B. Description of Slniulating Programs 
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The software for implementing the algorithms for Filter 3 and 
Filter 5 is developed at the assembly language level. The primary con- 
siderations in its development are the one byte word size and the lack 
of a hardware multiplier on the 6800 microprocessor. The consequences 
of the small word size are that internal values must be represented In 
multiple byte lengths. Since these values can only be written into or 
read from memory one byte at a time, execution, time is dependent on the 
number of bytes used to represent a value. For this reason, it is necessary 
to determine the minimum number of bytes needed to preserve precision at 
any given point in the program. Although the range input, range output 
and range rate output are two byte values, internal word Size varies from 
two to five bytes. 

Multiplications are handled by a multiplication subroutine. Under 
worst case conditions this routine requires 2,235 periods of the system 
clock to return a product. This represents a significant percentage 
of program execution time. 

1. Program FILTER 3 

An assembly language listing of the microprocessor 
realization of Filter 3 is presented in Appendix B. This program 
accepts a two byte unsigned binary word representing a noisy range 
value and outputs a two byte unsigned binary word representing 
the present estimated range and a two byte, two's complement 
binary word representing the present estimated range rate. 

When a noisy range value is received it is loaded onto the bottom 
of a 36 byte stack, REGS, containing the previous 17 noisy range 
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inputs « A summation is then made of the six most recent entries 
onto the stack and subtracted from that is the summation of the 
six oldest entries on the stack. The result of these operations 
is a three byte, two's complement number stored in the variable 

TSIIM. This value is multiplied by the scaled constant (1/(72*DT)) 

Ifi Ifi 

* 2 to determine the estimated range rate. The 2 scaling 

factor is included in the constant because the constant is typi- 
cally less than one. Since integer arithmetic is employed, 
tring to represent a number between zero and one by a zero or 
a one would cause a very significant loss of accuracy. To 
circumvent this situation the constant is premultiplied by a 
factor Of 2^®, which is equivalent to shifting the constant 
two bytes to the left. After the multiplication of TSIIM by 
the constant has been accomplished the two least significant 
bytes of the five byte product are truncated which is equivalent 
to a two byte shift back to the right. The magnitude of the 
estimated range rate is such that it can be represented suf- 
ficiently by a two byte value so the most significant byte of 
the three remaining bytes of the product can also be discarded. 

The end result is a two byte, two's complement number, ROOT, 
representing the estimated range rate, 

TSUM is multiplied fay the scaled eonstant (7,5/72)*2^® 
to derive a range correction termi, RGRCTj which will be used in 
the range ealGulations to give the estimated range. In this 
ease, the two lowest significant bytes are retained until after 
the addition in the range ealculations to preserve aecuraey. 
However, the most significant byte can still be discarded leaving 
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a four byte, two*s complement term RCRCT. 

To derive the estimated range, RE, the summation of the 
16 most recent entries onto the stack REGS is taken. The result 
is a three byte number located in SUM. This value is multiplied 
by the scaled constant (1/16)*2 producing a five byte result. 

The most significant byte is discarded and the remaining four 
byte value is added to RGRCT. The two most significant bytes of 
the sum from this addition are the unsigned binary representation 
of the estimated range. 

At this point FILTER 3 receives the next noisy range 
term and pushes it onto the bottom of the stack REGS causing 
the oldest entry on the stack to be Tost off the top. The process 
of obtaining the estimated range and range rate is then repreated 
for these 18 entries on the stack. 

Program FILTER 3 utilizes four subroutines: SHFTR, ADD, 

SUB and MULT. Subroutine SHFTR shifts all entries in the stack 
‘.EGS up one and loads the latest range error term onto the bottom 
of the stack. Subroutine ADD adds two three^byte numbers. The 
numbers are entered into the subroutine through the variables A 
and B. The three byte sum is returned through the variable SUM. 
Subroutine SUB is a three byte subtraction routine. Values are 
entered and returned through the same variables used in subroutine 
ADD. 

All multiplication required by FILTER 3 is performed by 
subroutine MULT, a two byte by three byte, two's eomplement 
multiplication routine. The two byte multiplicand is entered 
through MLTCD and the three byte multiplier is entered through 
MLTPR, The five byte product is returned through the variable 


PROD. 



INPUT TYPE OF PRN5E INPUT CRFiMP.- NUISE DR 
PHNP+NaiSE>. IMPUr LENGTH OF SINULRTIGN 
CNCI. OF S.fiMPL£5:> . SET UP POINTER FUR BLOCK 
OF PhNGE INPUT SRMPLES CRN> . SET UP POINTER: 
FOP BLOCKS OF OUTPUT SRNPLES CRRNGE-RE, 
iRPiNGS RRTE-RIiQT> , 


INPUT RM 


LOFtli SUN ONTO :STRCK ” REGS 


UNMFrriCN OF THE :SI:K NEUEST EHTIE:S If) RES 
NUS c:S:!JMMFfTlOr-i OF THE SIK □ 


TSUN ♦ <1 


RDOT c;N 


UN + < 


RDI CN.^ 


UNMfiTION OF The ip. NEUEST entries in REGS = 




OF OUiPU) 


IS SINJJLHrlDrT 
.^lOMPLETEIi? ^ 


STOP 


FIGURE 66 - FlQWGftart for program FILTER 3 
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Subroutine MULT also contains a rounding routine^ Since 

1 fi 

one of the factors of PROD has been scaled up by 2 , PROD must 

be scaled down by 2^® to obtain the correct result. Rather than 
truncating, PROD is rounded off to its three most significant 
bytes. The flag ROUND is used to turn the rounding routine on 
or off as needed. 

To facilitate the display of data the FILTER DATA PRINTER 
program is appended to the FILTER 3 program. For each sample, 
the printer program lists the sample number, the value of the 
noisy range input at that sample, the value of the actual range at 
that sample and the estimated range and range rate calculated by 
FILTER 3. This program converts all hexidecimal values to deci- 
mal values. 

2. Program FILTER 5 

An assembly language listing of the microprocessor reali- 
zation of Filter 5 is presented in Appendix G. This program 
accepts a two byte unsigned binary word representing a noisy 
range value. After receiving sixteen of these values the program 
outputs a two byte unsigned binary word representing the present 
estimated range and a two byte, two's compliment binary word 
representing the present estimated range rate. Also calculated 
but not output is a two byte unsigned binary word, RDI, represent- 
ing the estimated range at a previous point in time. This value 
is subtracted from each incoming noisy range word to produce a 
two byte, two's complanent binary rangie error term. These error 
terms are placed on a 32 byte stack, REGS. When sixteen noisy 
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range samples have been received, RIGS is then full and a sunwa- 
tion of Its contents Is taken. The result is a 3 byte; 
two's complement binary value, TSUM, which is used to calculate 
the new estimated range rate, ROOT, and to update the two estima- 
ted range values, Ri and RDI, 

The estimated range rate, ROOT, is derived by multiplying 
TSUM by the scaled constant (1/(256*DT))*2^®. Since the constant 
is typically less than one and integer arithmetic is employed, 
the scale factor 2 is included to maintain the two byte level 
of accuracy. The multiplication of the 2 byte constant and 3 
byte TSUM results in a 5 byte product. The estimated range rate 
value is contained in the three most significant bytes of the 
product due to the scale factor included in the constant. This 
term can be accurately expressed in two bytes so that the most 
significant byte can also be discarded yielding a two byte, 
two's complement, estimated range rate value. 

The three byte TSUM is then multiplied by the scaled 
constant {l/i6)*2'*^ and the resulting five byte product is 
reduced to two bytes in the same manner as described above. 

This two's complement number is added to the present value of 
RDI to update it for use in determining the range error terms 
for the next sixteen noisy range inputs. RDI is a two byte, 
unsigned binary value. 

The estimated range, RI, is updated In the sane manner. 

TSUM is muliptlied by the scaled constant (23.5/256) *2^®. The 
resulting five byte product is reduced to a two byte, two's 
complement number and added to the old value of RDI, The result 
of this addition is a two byte, unsigned binary number representing 
the current estimated range. 


INPUT TYPE OF RANGE INPUT <F:RNP!. NOISE OR 
RRMP+NOISE>. INPUT LENGTH OF SIMULPTIPN 
CNO, OF SRNPLES? . SET UP POINTER FOR JE:LDCK 
OF RANGE INPUT SAMPLES *;;RM> . SET UP Pu INTER;; 
FOR BLOCKS OF OUTPUT SAMPLES fiRANGE-RE.i 
RANGE RATE-RriOT> . 


INPUT RM 


RM - RUI = SUM 


LOAD SUM ONTO :STACK RE( 


HAVE 16 INPUT 
BEEN RECEIVIIi 


SUMMATION OF CONTENTS OF RE( 


TSUM C 


6 ^ It:':- = RJitOTCN.' 


asuM 


+■ RBiaC- = RHI CN+1> 


<TSUM ♦ «;E3.5.'-£56>:> + RHI CN> = RECN:-' 


is STMULATIOfT 
OCOMPLETEB? ^ 


STOP 


FIGURE 67 ^ Flowchart for program FILTER 5 
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Program FILTER 5 utilizes the same four subroutines used 
by Filter 3» The operation of these subroutines is exactly as 
described in the section on FILTER 3. Temporary storage is allo- 
cated for the stack REGS in such a way that subroutine SHFTR, 
which is responsible for pushing data onto the stack and shifting 
the other entries up, need not be changed to reflect the fact 
that REGS is a 36 byte stack in FILTER 3 and a 32 byte stack in 

FILTER 5. Both programs utilize premul tiplication of donstants 
1 

by 2 so that the rounding routine in subroutine MULT also need 
not be changed. 

Since FILTER 5 outputs estimated range and range rate 
data only once every sixteen sample periods, the FILTER DATA 
PRINTER program is modified to display data for only those sample 
periods for which outputs are generated, 

G, Test Conditions 

As noted previously, test conditions for the microprocessor 
simulation runs are very similar to those used in the IBM 370 simula- 
tions. Tests are conducted for each of the three range input signals; 
the ramp, Gaussian noise and the summation of the ramp and noise. The 
tests are hot conducted in real-time but rather a block of data eontainihg 
sequential samples of the input is resident in memory. The simuTatihg 
program steps through this block retrieving input values as needed. A 
block consists of 630 sampled range values taken at the sampling rate. 
Each value is a two byte uhstgned binary word. These are rounded integer 
representationsof the same sampled input values employed in the IBM 370 
simulations. 
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The microprocessor simulations are conducted at the 6.25 samples/ 
second sampling rate only. Although simulation at the higher sampling 
rate is possible, the speed limitations of the microprocessor would make 
a real-time system unrealizable^ The two filter designs chosen both 
offer better noise rejection at the lower sampling rate. 


D. Test Results 


1, Filter 3 

Simulation runs for Filter 3 are conducted for a duration 
of 200 sample periods. The last 50 sample values of the range 
input, estimated range and estimated range rate are statistically 
analyzed and presented in Tables XLIX - LI. Also presented is 
the data for an equivalent run on the IBM 370. The results indicate 
that the implementation of the filter on a microprocessor based 
system is accomplished without significant loss of accuracy. 

Program FILTER 3 requires 15,777 cycles of the micropro- 
cessor system clock to compute estimated range and range rate 
values after being given a sampled range input. At the 614.4KHz 
clock rate of the MEK 6800D2 Evaluation Kit, the execution time 
is 25.68 milliseconds. The 160 millisecond sampling period provides 
more than adequate time for filter Galeulations, hence, this program 
could be employed in a real-time system. 


2. Filter 5 

Filter 5 simulations are cOndueted for a duration of 
608 sample periods. Since this filter returns an output only once 



^Hn mm mm 


■ i 


TABLE XL IX 

FILTER 3 RANGE STANDARD DEVIATIONS(FT, ) AT 6.25 SAMPLES/SEC. 


FILTER 3 

RUM 1 

input = ramp 

RUN 2 

input = noise 

RUN 3 

imput = ramp + noise 

INPUT 
RANGE 
SliND. DEV. 

OUTPUT 

RANGE 

STND;. DEV. 

INPUT 

RANGE 

STND. DEV. 

* 

OUTPUT 

RANGE 

iSTND. DEV. 

INPUT 
RANGE 
STND. DEV. 

OUTPUT 
RANGE 
STND, DEV. 

M1GROPROCESS0R 

^ 65.380 

65.375 

103.531 

53.105 

102.647 

61.596 

IBM 370 

65.3068 

66.3068 

105.9378 

53.1095 

108.1119 

68.1752 

■ FILT! 

FILTER 3 

ER 3 RANGE RATE STANDARD 
RUN 1 

iinput * ramp 
OUTPUT 
RANGE RATE 
STND. DEV. 

TABLE L 

0EVIATI0NS(FT./SEG.) AT 6. 
RUN 2 

i■n;^ t = noise 
5TPUT 

. RANGE RATE 
STND. DEV. 

25 SAMPLES/SE 
RU 

input = ra 
OUT 
RANGE 
STND. 

:c. 

H 3 

rnp + noise 
PUT 
RATE 
DEV. 

MICROPROCESSOR 

0.0 

33.288 

33.132 

IBM 370 

0.0 

33.1917 

33.1917 














TABLE LI 


FILTER 3 RAHGE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEC. 


FI LTER 3 



2. 632 X 

1.034 X 10“^ 

3 

2.513 X 10"^ . 

9.817 X 1'0~^ 

3 
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for every sixteen sampled inputs, 38 output values are deter- 
mined during a simulation run. The filter converges to steady 
state after two execution cycles (32 sample periods), therefore, 
the last 35 output values represent steady state and are statistically 
analyzed. This set of samples is of sufficient size to insure 
statistically accurate results. Data for the microprocessor 
simulations of Filter 5 and an equivalent IBM 370 run are presented 
in Tables LIT - LIV= Here again, as in the microprocessor reali- 
zation of Filter 3, results agree with those obtained on the 
IBM 370, 

The worst case execution time for an "OFF" sample period, 
when no filter outputs are produced, is 907 cycles of the micro- 
processor system clock. The execution time for an "ON" sample 
period, when filter output values are determined, is 11,639 
cycles of the system clock. At the 614.4KHz clock rate of the 
MEK 6800D2 Evaluation Kit, execution times are 1.48 and 18,94 
milliseconds for "OFF" and "ON" sample periods, respectively. 

However, the 6800 family is capable of clock rates up to 2MHz 
which would yield execution times of 0,45 and 5.82 milliseconds 
for "OFF" and "ON" sample periods, respectively. 

To offset the problem of smoothed range and range rate 
data being available only every sixteenth sample instant, it 
has been proposed that sixteen filters be operated in parallel 
with each filter lagging the next by one sample period. In 
this manner an output from one of the filters would be available 
at each sample instant. At the 6.25 samples/second sampling rate 









TABLE LM 

FILTER 5 RA'NGE STANDARD DEVIATIONS{FT. ) AT 6.25 SAMPLES/SEG. 


FliLTER 5 

RON 1 

IriDUit = ramo 

RUN 2 

in out = noise 

RUN 

Inouit = ram 

3 

+ noise 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV, 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

STND. DEV. ^ 

INPUT 

RANGE 

STND. DEV. 

OUTPUT 

RANGE 

^ STND. DEV. 

microprocessor 

734.514 

734,516 

91.516 

32.489 

755.846 

735.662 

IBM 370 

734.5007 

734.5007 

111.0584 

35.4184 

739.8486 

735.4358 

FILT 

i 

FILTER 5 

ER 5 RANGE RATE STANDARD D 
RUN 1 

input = ramp 
OisITPUlf 
RANGE RATE 
STN13. DEV. 

TABLE LI II 

EVIATIONS(FT./SEC.) AT 6.2 
RUN 2 

Input = noise 
OUTPUT 
RANGE RATE 
STND. DEV. 

5 SAMPLES/SEC. 

RUN 3 

input = ramp + noise 
OUTPUT 
RANGE RATE 
STND. DEV . 

MIGROPR0GESS0R 

0.0 

11.543 

11.508 

HBM 370 

0.0 

12.5710 

12.5710 
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TABLE LIV 

FILTER 5 RANvSE AND RANGE RATE VARIANCE REDUCTION RATIOS AT 6.25 SAMPLES/SEG 
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the execution time of program FILTER 5 is only a fraction of the 
sampling period. Hence, the same processor can be used for all 
sixteen filters. In a given sampling period fifteen filters 
would be at “OFF" sampling periods while one would be at its 
"ON" sampling period. Operating with the 614,4KHz system clock, 
total execution time would be 42.62 milliseconds, slightly over 
a quarter of the sampling period. 


CONCLUSIONS AND RECOMMENDATIONS 


The IBM 370 simulations demonstrate that Filter I offers im- 
proved noise reduetion at higher sampling rates in both the range and 
range rate modes for both constant gains and band! imi ted differentiator. 
Maneuver- following capabilities, as measured by convergence time, are 
independent of sampling rate. Noise reduction for Filter 2, Filter 3 
and Filter 5 is independent of sampling rate in the range mode and in- 
creases with decreasing sampling rate in the range rate mode. These 
filters all demonstrate increased maneuver-following capability at 
higher sampling rates. Filter 4 does not provide sufficient noise re- 
duction to be considered further. 

At the low sampling rate (6.25 samples/second) and with a short 
Gonvergence time. Filter 5 has greater noise reduction than Filter 3. 

If parameter pairs are Ghoseri for Filter 1 and Filter 2 to yield an 
equally short convergence time, they too do not provide the level of 
noise reduction of Filter 5. However, if a longer convergence time is 
allowed, the parameter pairs of Filter 1 and Filter 2 can be chosen to 
provide higher levels of noise reduction than Filter 5. In this case. 
Filter 1 (bandltmited differentiator) provides greater noise reduction 
than Filter 2 for a given convergence time. 

The advantage of Filter 1 becomes more apparent at the higher 
sampling rate. Convergence time is unaffected but noise reduction 
increases with sampling rate. Hence, the caipablTities of this filter 
are determined by the- limitations of the processor with which it is 
implemented. Future generations of proeessors will offer faster clock 
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rates which will allow a corresponding increase in sampling rates and 
result in better filter performance. 

Filter performance can be summarized: 

Filter 1: t - independe*'t of sampling rate 

- dependent on Kl, K2 

KR - decreases with increasing sampling rate 

- dependent on Kl, K2 

KRR - decreases with increasing sampling rate 

- dependent on Kl, K2 

The tradeoff between t and the ratios KR and KRR is controlled 
by Kl, K2 (b in the BLD). 

Filter 2: t - decreares with increasing sampling rate 

- dependent on a,B 

KR - independent of sampling rate 

- dependent on a,g 

KRR - decreases with decreasing sampling rate 

- dependent on a,B 

The tradeoff between t and the ratios KR and KRR is controlled 
by cx,Bx 

Filter 3: t - decreases with increasing sampling rate 

KR - dependent on sampling rate slightly, due to updating 
KRR - decreases with decreasing sampling rate 

Filter 5: x - decreases with increasing sampling rate 
KR - independent of sampling rate 
KRR ^ decreases with decreasing sampling rate. 

Future study should be directed toward implementation of Filter 

1 or Filter 2 at a higher sampling rate with a 16 bit mieroproeessor 

having a hardware multiplier. It would also be worthwhile to investigate 

different methods of approximating the continuous system of Filter 1 

with discrete systems. The method used in this study is only one way of 

mapping from the s-plane to the z-^plane. Other methods should be examined 

to determine if any offer mappings more suitable for this application. 
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APPENDIX A 


The follGwing pages contain a FORTRAN listing of program FLT, 
the associated GALG subroutines for each of the five filters and the 
two plotting subroutines, PLOTS and MGPL0T= 


o n o 



c 

r THIS program S:IMIJ!LATES a digital filter. 

C THE PARAMETERS OF THE FILTER AND IMITIAL CONDITIONS ARE READ IN 
C ON NAMELIST LIST. 

C THE PLOTTING PARAMETERS ARE READ IN OM NAMELIST OTM. 

C THE PARAMETERS OF THE RANDOM NUMBER GENERATOR ARE READ IN CN 

c naMielist xrand. 

C SUBROUTINE CALC PERFORMS ALL FILTER GALCULATIONS. 

C SiU'B ROD TINE GGNRF GENERATES THE NORM.ALLY DISTRIBUTED RANDOM VARIABLES 
C THAT ARE USED AS RANOE NOISE. 

C SUBROUTINE BECORl eALCULATES THE MEAN AND STANDARD DEVIATION OF THE 
C INPUT RANiGE* output RANGE AND OUTPUT RANGE RATE. 

C SUBROUirNE PL 0TB FORMATS DATA TO USE THE PRINTER FOR PLOTTING. 

G SUSROOTJNE HOP LOT FORMATS DATA FOR THE GALCO'MP PLOTTER. 

C 

C PROGRAM FLT MAKES THREE FILTER SIHOLATION RUNS FOR EACH SET OF 
C FILTER PARAMETERS. A NEW SET OF FILTER PARAMETERS IS READ IN EVERY 
C TIME NAMELIST LIST IS READ. THE FIRST SIMULATION RUN IS MADE WITH A 
G RAMP INPUT. DATA FROM THIS RUN IS TAB'ULA TED IN A TABLE ONLY (NO 

C PLOTS). the SECOND SIMULATION RUN IS MADE WITH AN INPUT OF GAUSSIAN 

C NOISE. DATA FROM THIS RUN IS ALSO TABULATED IN A TABLE ONLY. 

G THE THIRD SIMiULATION RUN IS MADE WITH THE INPUT OF A NOISV RAMP, 

e DATA FROM THIS RUN IS TABULATED IN A TABLF AND ALSO PLOTTED. 

C THE PLOTS RETURNED ARE DETERMINED BY THE PARAMETERS IN NAMELIST DIM. 

G ALL PLOTS ARE OF A FUNCTION VERSUS TIME. 

O 

G NAMELIST LIST 

G LIST IS READ AT THE BEGINNING OF EACH RUN THROUGH THE MAIN PROGRAM 

t DiRE = INITIAL ESTIMATE OF RANGE 

C DRDOT = INITIAL ESTIMATE OF RANfGE RATE 

DT = TIME BETWEEN INPUT SAMPLES 
K1 AND K2 = PARAMETERS OF FILTER 
DR ,= INITIAL RANGE 
C DV = INITIAL RANGE RATE 

-p. 
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MDM = NUMBER OF FILTERS TO B£ SIMULATED <THE NUMBER OF TIMES 
NAMELIST LIST IS TO BE READ.) 

IDl = NUMBER OF SAMPLE PERIODS BETWEEN THE READING OUT OF 
RESULTS 

IDRI = NUMiBiER OF RESULTS TO BE READ OUT 

NAMELIST DIM 

DIM IS READ OiNLY OWCE. 

PER = 1 SPECIFIES THE OUTP'UT OF TWO RANGE PLOTS, BOTH 

PLOTS CONTAIN THREE FUNCTIONS-R »R M AND RE, THE FIRST 
PLOT USES A SCALE CHOSEN FR CM THE MAX AND MIN FU'NGTION 
VALUES. THE SiECONiD USES BUI AND BLl AS THE PLOT'S MAX 
AND MIN VALUES. IF PLR IS NOT SPEGiFIED NO RANGE PLOTS 
ARE GIVEN, THE PLOTS ARE GENERATEO AFTER THE THIRD 
FILTER SIMULATION RUN FOR EACH FILTER, 

PLRR = 1 SPECIFIES THE OUTPUT OF TWO RANGE RATE PLOTS, BOTH 
PLOTS CONTAIN TWO F#NGTIONS-V AND RDOT, THE FIRST 
PLOT USES A SCALE GHOSEW FROM THE MAX AND MIN FUNCTION 
VALUES, THE SECQ'ND USES B U2 AND BL2 AS THE PLOT'S MAX 
AND MIN VALUES, IF PLR IS NOT SPEGIFIED NO TANGE RATE 
PLOTS ARE GIVEN. THE PLOTS ARE GENERATED AFTER THE 
THIRD FILTER SIMULATION RUN FOR EACH FILTER. 

PLRER - 1 SPECIFIES THE OUTPUT OF ONE OR TWO RANGE ERROR 

PLOTS. UNLIKE THE FIRST TUG SETS OF PLOTS* THESE ARE 
NOT GENERATED UNTILL ALL FILTERS WAVE BEEN TUN. 

HENCE* THE RANGE ERROR FOR THE THIRD SIMULATION RUiN 
FOR EACH FILTER APPEARS ON THE SAME PLOT. THIS IS 
FOR COMPARISON PURPOSES. THE FIRST 

PLOT USES A SCALE CHOSEN FROM THE FUNCTION'S MAX AND 
MIN VALUES, THE SECOND USES BU3 AND BL3 AS THE 
PLOT'S MAX AN® MIN VALUES. IF PLRER IS NOT SPECIFIED 
NO R A N'6 E E R ROT PL 0 T S AR E 6E N£ R A TE D . 

PLRRER = 1 SiPECIFIES THE OUtPUT OF ONE OR TWO RANGE RATE 

ERROR PLOTS, THESE ARE ALSO NOT GENERATED UNTIL ALL 
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FILTERS HAVE RFEW RUN. HENCE, THE RANGE RATE ERRCP 
FOR THE THIRD SIMULATION RUN FOR EACH FILTER APPEARS 
ON THE SAME PLOT. THE FIRST 

PLOT USES A SCALE CHOSEN FROM THE fUNCTION*S MAX AND 
MIN VALUES. THE SECDWD USES BUA AN0 BL4 AS THE 
PLOT’S MAX AND MIN VALUES. IF PLRRER IS NOT 
SPECIFIED NO RANGE RATE ERROR PLOTS ARE GENERATED. 

BUI AN® BLl = THE UPPER AM® LOVjER LIMITS ON THE SECOND GRAPH 
OF THE RANGE DATA. If NOT SPECIF lED THEY ARE 
SET 

BUI = (MEDIUiM ACTUAL RANGE) ♦ SPAN! 

B'Ll = (MEDIUM ACTUAL RANGE) - SPAM 

U'MCRE SPAN! IS A MAIN PROGRAM GONSTANT. 

BU2 AND BL2 = TMiE UPPER AND LOWER LIMITS ON THE SECOND GRAPH 
OF the range rate DATA. IF NOT SPECIFIED THEY 
ARE INTERNALLY SET 

BU2 = (MEDIUM ACTUAL RANGE) + SPAN2 

BL2 = (MEDIUM actual RANGE) - SPAN2 

WUERE SPAN? IS A MAIN PRCGRAM CONSTANT. 

BU3 AND BL3 = THE UPPER AN® LOWER LIMITS ON THE SECOND PLOT 

OF THE RANGE ERROR FUNCTION. IF NOT SPECIFIED, 
H SECOND PLOT IS NOT GENERATED. 

BU4 AND BLA = THE UPPER AND LOWER LIMITS ON THE SECOND PLOT 
OF THE RANGE RATE ERROR FUNCTION. IF NOT 
SPECIFIED, A SECOND PLOT IS NOT GENERATED. 

GALGOM = 1 CREATES CALCOMP PLOTS. A RANGE PLOT AMD RANGE 

RATE PLOT AFTER THE THIRD SIMULATION RUN FOR EACH 
FILTER. THE RANGE ERROR AND RANGE RATE ERROR 
PLOTS ARE generated AFTER ALL FILTERS HAVE BEEN 
RUN, HENCE, THESE PLOTS CONTAIN THAT DATA FROM 
THE THIRD SIMULATION RUN FOR EACH FILTER. THIS 
IS FOR COMPARISON PURPGSES. 
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C *N0TE* INTERMA;L SPECIF IGA TIGN 0F BUl»8LlrB'U2 AND BL2 CAN 

G G'NLY RE USED WHEN NO AGGELERATION IS A&DED. IF THE 

G VALUES OF SPANl AND SPAN2 ARE NOT ADEQUATE, THEY HAY 

'BE GHANiGED BY ALTERING THEIR VALUE IN THE MAIN 
PROGRAM* 

NAMEsLIST XRANI0 

XRANID IS READ ONLY ONCE • 

XM 3 THE MEAN OF THE NORMALLY DISTRIBUTED RANDOM VARIABLES 
XDEV = THE STANDARD DEVIATION O'F THE NORMALLY DISTRIB'UTEO 
RANOO'M VARIABLES 

INTERNAL PARAMETERS 

NiN = G TO NDN, NUMiBER OF FILTERS THAT HAVE BEEN SI HUE A TED. 

L = 0 TO IDRI , NUMBER OF TIMES SURROUTINE CALC IS CALLED 
e (ONE DATA POINT IS GENERATED FOR EACH L> 

C 

G 

C PROGRAM FLT 
C DOUBLE PRECISION 

IMPLICIT REAL*BtD,K,R,V) 

I N T E G€ R P L R , P L R R , P L RE R , F L RR E R , R U'M 
DIMENSION TIMA (500) 

DIMENSION X ( 50 ,5 ),DIFF( 5,500) ,RERR (5, 500 ) ,Y (2500) 

D I ME NS I O'N R NfG ( 5 , 5 0 0 ) ,R NG R ( 5 ,50 0 ) , T I ME ( 500 ) 

DIMENSION SNG(5,50O) ,SNGR(5,500) ,SIFF( 5,500) ,SERR(5,500) 

D I ME N S I OIN 2 ( 6 , 6 ) , X M I ( 6 ) , S D I < 6 ) 

G MATIX 7 . IS A DiUMiMY VARIABiLE USED IN BE COR I. 

N A ME L I S T V L I S T /OiR E , 0 R D O T , D T , K 1 , K 2 , DR ♦ D V , N® N » I D I , I D R I 
N,A M E L I S T / D I M /P'LR , P LRR , P L RE R , P LR R E R , 0 U 1 ,B Ll , BU 2 , B L2 , B U 3 , 31 3 , B UA , B L4 
2, GAL COM 

WAMCLIST /XRANO/XH ,XDEV 
CALL PLOTS (lG2A,<>a96) 


VO 


CALL PtTHSG (•*** BON FLEISCHER PLOTS READY ***««33) 
CALL PLOT (0.0*- 11. 0,-3) 

CALL PLOT (0.0,1.0,-3> 

W'N = 0 

BtJl =0.0 

BLl =0.0 

BU2 = 0.0 

BL2 = 0.0 

SU3 = 0.0 

BL3 = 0 .0 

tUA = 0.0 

BLA =0.0 

PLR = 0 

PLRS = 0 

PLRER = 0 

PLRRER = 0 

CALGOM = 0.0 

GALL GGNRFl (DSEEB> 

READ (5, DIM) 

READ(5fXRANID) 

100 READ(5,LIST) 

MM = NN+1 
BiETA = I0I*DT 
110 RUN = 1 
RE = ORE 
ROOT = DlRDOT 
R = DR 
V = DV 
XXM = XM 
XXOEV = XOEV 
X M = 0 . 0 
XDEV = 0.0 
GOTO 1A0 
RUN = 2 


12 0 
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CALL GgN'Rf 2(0 SEED) 

RE = 0.0 
ROOT = 0.0 
R = 0.0 

V = 0.0 
X« = XXH 
XOEV i XXiDCV 
GOTO 140 

130 RUN = 3 

GALL GGN‘RF2< DiSEED) 

RE = 0iRE 
ROOT = DR DOT 
R = DR 

V = Dtf 

140 eONTINUlE 

WRITE(6*150) 

150 FORMAT C ItHOf 2X* 'LIST* ) 

W R I T E ( 6 » I fe O ) R E * R D 0 T ♦ D T , K I , K 2 t R , V * N'D N ♦ I D I , I DR I 
160 F0a»AT(lH0r4X,*RE = • t G2 0 .5 /5X , * R DOT = • ,G20 .9/5X , »0T = **G20.9/5X 
2»»K1 = * ,G2D.9/5X**K2 = • , G2 0 ,R/5X ,* R = • ,G20 .9/5X , » V = •»G20.9/5X 
3»»NDN = •,I11/5X,*!DI = NIll/SXrMDRI = NUl///) 

WRITE(fi*170) 

170 FORHAT(lHO,2Xt*XRAN!D«) 

WRITE(G*ieO) X'M,XDEV 

160 FQRMATUN0*4X,*XM = * ,G20.9/5X , *XDEV = •»620.9) 

IF (RUN .ME. 3) GOTO 300 

GALeULATION OF B#1.8Ll,iBU2 ANiD BL2 IF NOT GIVEN AS A INPUT 

IF (BUI .N£. 0,0 .OR. BLl ,NE. 0.0) GOTO 190 
SPANl = lOOO.G 

BUI = R ♦ V*I:D'I*DT*I0RI/2 ♦ SPANl 
BLl = BUI - 2*SPAN1 

19Q IF (BU2 .NE. 0.0 .OR. BL2 .NE. 0.0) GOTO 200 
SPAN2 = 2.0' 


TGT 
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iU2 = V + SPAN2 
BL2 = V - SPA^'2 
200 WRITc;(6,2lD> 

210 F OiR M A T n H 0 , 2 X * • D I M • ) 

URITE<6 ,2201( PLR »PLRR tPlRER iPLRRERtBUl »BL1»BU2»BL2 
220 FOiRMAT(lMO»ftX,*PLR = • , 1 6 /5 X , • PLR R = * » I 6/5 X , • P LRER = »tl6/5X,*PLR 
2RER = • ,I6/'5X»*BU1 = • » G2 0 . 9/5 X , • B LI = * , G2 D ,9/5X , »BU2 = *<620 . 9/5 
3X»*BiL2 = *tG20.9) 

IF (RiU3 .EQ, O.G .AND, BL3 .ID. 0.0) GO TO 240 
WRITE (6, 230) BU3*BL3 

230 F0RMiAT(lH0,4X»*BU3 = ♦ ,62 0 .9/5X » *BL3 = »,G20,9) 

GOTO 260 

240 WRITE (6,250) 

250 F0RMAT(1H0,4X, *BU3 AND 8L3 WERE NOT SPECIFIED IN DIM») 

260 IF (BD4 iEQ. 0.0 .AND. BL4 .EQ. 0.0) GO TO 280 
WRITE (6, 270) B.U4,Bi4 

270 FORMAT! 3H0,4X,«BD4 = • , G2 D . 9 /5X , * BL4 = *,G20.9///) 

GOTO 300 

260 WRITE (6,290) 

290 FOWMAT(1HO,4X ,«BiU4 AND BL4 WERE NOT SPECIFIED IN DIM*///) 

300 CONTINUE 
L = 0 

WRITE! 6,310) 

310 F OiRM AT ( 1 HO , 2 X , • T I ME • , 1 7X , * R • , 1 7 X , * RE » , 15X , * R ERR • , 16 X , • V ♦ , 1 5X , • R D OT 
2* ,15X,*DIFF* ) 

RDI = 0.0 
RVI ~ 0.0 

e 

G AN IF STATEMENT WILL BE INSERTED HERE TO CHECK TO SEE IF THE RANDOM 
C NUMBER GENERATOR IS TO BE USED ON THIS RUN. IF IT IS, THE FOLLOWING 
e LOOP WILL BE SKIPPED. 

C 

320 IF (XDEtf .ME. 0.0) GO TO 340 
DO 330 LAP = 1,IDI 
Y(LAP) = 0.0 
COWTIiMiUE 


330 


tn 


GOTO 350 


I 


C 

C THE GALL STATEMENT FOR THE RANDOM N, UMBER GENERATOR WILL RE 
C HIEKE. THIS SUBRCUTINiE WILL BE CALLED EACH TIME CALC IS CALLED. 

e 

340 GALL G'GNRF nDI«Y) 

DO 350 LAP = 1*IDI 
V(LAP) = XM + X0EV*Y(LAP) 

350 eONTINUE 
L = L + 1 

TIHEILJ = IDI*DT*L 

CALL GALG(R,V,OT^X,IDI»RE<ROCT»K1<K2,R0I»RVI,L,RH,Y, IDRI) 

C 

C VARIABLES BEGINNING WITH R OR D ARE DOUBLE PRECISIGN, THEY ARE USED 
G IN DATA TABLE OUTPUT. VARIABLES BEGINNING WITH S ARE SINGLE 
C PRECISION. THEY ARE USED IN THE PLOT OUTPUT. 

C 

C R'NG GOMTAINS 3 FUNCTIONS TO BE PLOTTED ON THE RANGE GRAPH. 

RNGU»L> = RE 
RNG(2»L) = R 
BNG(3,L> = RM 
SNG<1,L) = RE 
SNG(2,L) = R 
SNG(3tL) = RM 

C R’MGR CONTAINS 2 FUNCTIONS TO BE PLOTTED ON THE RANGE RATE GRAPH. 
RNGRUrL) = ROOT 
RNGR(2»L) = V 
SNGR(lvL) = ROOT 
SN6R(2rL) = V 

e RERR CONTAINS 1 FUNCTION TO 3E PLOTTED ON THE RANGE ERROR GRAPH. 
RERR(NN»L> = RE - R 
SERRINNrD = RE - B 

G OIFF CONTAINS 1 FUNCTION TO BE PLOTTED ON THE RANGE 
C RATE ERROR GRAPH. 

DIFFCNNrL) = ROOT - V 


cn 

CJ 


SrFF<WN»L) = RDOT - V 
C 

W'R I T€ ( 6 « 3 60) T I MC ( L ) i R t R E ♦ R E R R ( M N * L > i V ,R D O T , D I F F ( N N , L ) 

3 6'0 F OR MA T ( 1 HO * 3 X *G7 • 1 »6 1 1 X , G 1 0 .11)) 

IFfi.LT.lORI) 60 TO 320 
R ERR MX = 0.0 
0IFFHX =0.0 
00 370 LAP=1,50 

IF{A8S(X(LAP*<H) .GT .RERRMX) RERRMX = A BS (X (L AP »4 ) > 

IF( A.BS(X(LAP*5) ) .6T .DIFFHX ) 0TFFMX = A ES ( X ( L AP , 5) ) 

370 CONtlMU'E 

W*RIT£( 6i3R0> 'RERRMX tO I FFMX 

380 FORM'ATtlHD»RX,*MAXIMUM STEADY STATE RAAGE ERROR = * »62 0 .9/5X , • MAX I 
inUM STEAiDY STATE RAMiGE RATE ERROR = '*G20.9) 

CALL BECORICX»50*3*50»XMI,SDI »Z,IER) 
m I T€ < 6 » 39 0) XM I I 1 ) .SO I ( 1 ) 

390 FORMAT(lHO,AXf»ir.)PUT RA^JGE MEAN = • ♦ G2 0 . 9/5 X , • I NPUT RANGE STND DEV 
1 = •,G20.9) 

yiRlTE (6.AOO) XMK2) .SDK 2) 

400 FORMAT <1 HO, 4 X* •range MEAN = • ,G2 0.9/5X ,* RANGE STND DEV = »,G20.9) 
WRITE (6, 4 10) XMK3) ,S0! <3) 

410 FDRHAT<1H0,4X, ‘RANGE RATE MEAN = •, G20 .9/5X , ‘RANGE RATE STND DEV = 


‘,620 

.9) 




IF(RUN 

•E 0 . 

1) 

GOTO 

120 

IF(R:UIN 

• EQ. 

2) 

GOTO 

13 0 

IF(PLR 

.EQ. 

0) 

GOTO 

460 


CALL PL G TB ( 3 , I DiR I , BE T A , S NG , T I M E , 0 . 0 , 0 . 0 , 0 ) 

•iRITE(6,420) 

W:RITE<6»430) 

WRITE (6,440) 

UR I TEC 6, 450) 

4 20 FO.RMAT(1HO,2 6X, *****, ‘FOR RANGE DATA* , 24X , ‘ * * * • ) 

430 F0)RMAT(1HO,2GX,‘***‘,‘1 IS THE CALCULATED RANGE ( R E ) ‘ ,S X, ‘ ** * • ) 
440 FORHATC IM0,26X, “***,»2 IS THE ACTUAL RANGE ( R ) * , 1 3X , ‘ ♦** • ) 

450 FGRMAT(1HD.26X,‘****,‘3 IS THE NOISY RANGE ( R M )*, 1 4X ,****• ) 








m 


^ »» 


cn 




CALL PLOTR ( 3 , 1 OP I , RETA , SNG , TI Mt ♦ R'Ul , B L1 1 0 ) 

URlT€(fe»420) 

L'^RITE(6»43tl) 

WRITE! 6»440) 

WR ITE! It, ,450) 

46fl IF(PLRR 4£Q. 0) GOTO 50G 

CALL PLOTS { 2,IDRif BETA,SNGR,TIMC,0.0, O.fO) 

URITE<6t470) 

WRITE (6 4480) 

WRITE (6 *490) 

470 FORMAT! IH0,26X,*****»*FCR RANGE RATE 0 AT A • » 1 9X * • * * * * ) 

4B0 F0RMAT(1R0,2F,X» •****f*l IS THE CALCULATES RANGE RATE (ROOT) * « 1 X , • * 
2 ** » ) 

490 FORMAT!lHO»26Xi*******2 IS THE ACTUAL RANGE R ATE ' ( V ) • ♦£ X , • * t > 

CALL PLQTB ( 2 ♦ I DR I , BET A t S NGR » T T ME t BU2 t BL2 1 0 > 

WRITE (6, 470 ) 

WRITE (6 ♦48 0) 

WR ITE t 6t490) 

50 0 CO NT I NIUE 

IF(CALCOM .EQ. 0.0) GOTO 510 

CALL f^l C PL DT ( S NG « T I ME ♦ I D R I ♦ 2 ♦ 0 ♦ $ W , D A T A ) 

GALL MCPLDTlSN'GiR .TlMEt ID'Rl»2»0fSW»DAT A ) 

I DR 1 1 = IDRI - 15 
DO 50 5 I = l.rI0iRII 
I A = I + 15 
TIMA (I) = TIMET I A) 

SNG(l.I) = SNG(1»IA> 

SN,G{2*I) = SNG(2,IA) 

SNGR{l,I) = SNiGR(lflA) 

SNiGR(2fI) = SNGR(2,IA) 

505 CONTINUE 

GALL MePLOTtSNG»TIMA»IDRII%2>0,SU*DATA ) 

CALL HGPLOTT SN'GR ♦TIMA,ID'RI1 ♦2,0,SW*DATA) 

510 IF(NN .LT. NDW) GOTO 100 

IF (PLRER .EO. 0) GO TO 530 
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CALL PLOTB ( WN * I DR I * B£ T A , SER R » T I KE , 0 . 0, 0 ,0 , 1) 

WiRITE t& ,5 20) 

52D FGRMAK lHO,2f,X,*****,»fOR RANGE ERROR LR ERR }», 1 GX ,****• ) 

IF (BU3 ,E9. 0.0 .AND. aL3 .EQ. 0.0) GO TO 530 
CALL PL 0TB t NiDN 1 1 DiR I * BiE T A t S ERR j T I HE , B U 3 * BL 3,1) 

WRITE ( 6, 520) 

530 IF (PLRRER .EQ. 0) GO TO 550 

CALL PLCTR ( NiDN , J D'R I , BiETA ♦ S I FF ,T I ME, 0 . 0,0 . 0 , 1 ) 

WiRITE(6,5ftO) 

5A0 FQiRMAT(lMQ»26X, ****», 'FOR RANiGE RATE ERROiR ( DI FF) • ,1 IX ,****• ) 
IF (BU'A .EQ. O.Q .ANiD. BIA .EQ. 0.0) GO TO 550 
G A L L P L 0 TB ( NID'N , I D R I ♦ B E T A ♦ S I F F , T I HE , B U A, B L A , I ) 

WRITE (6, 5 AO) 

55 0 COiNTTNiUE 

IF (GAL COM .EQ. 0.0) GOTO 57 0 

G THIS LGQP GENERATES A ZERO REFERENCE LINE FOR THE CALCOpp PLOT. 
NiDNl = MDN ♦ 1 
DO 5fe0 L = 1,1 DR I 
SERR(NDN1,L) =0.0 
560 SIFF (NiDNl,L) = 0.0 

e A LL HIG P LO T ( SE R R , T I ME , I D R I , NDN I , 0 tSU ♦ D AT A ) 

GALL HGPLGiT(SlFF,7IME,IDRI,NDNl ,0,S«,DATA) 

DG 56 7 0.= 1, N0IN1 
DO 5 65 1=1, I DR II 
lA =1+15 

SIFFtJ,!) = SIFF(J,IA) 

SERR(d, I) = SERR(d,IA) 

5 65 CO NT I NUE 
567 CONTINUE 

CALL MCPLGT (SERRrTIM'A, IDiRlI,NiD'Nl,0,SW, DATA) 

CALL M'CPLOK SIFF,TIMA,ID:RII,NiDNl,0,SW,DATA) 

570 CD NT I NUE 

GALL PLOT (0.0,0.0,999) 

STOP 

ENO 


9GT 
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SUf ROliT INE P LGTE { NP , MA X ,BO A « FNCT ,T IME, BU, BL « NZERO ) 

, PLOTB USES TWE PRINTER AS A PLOTTER. 

NP = KUMBER or FUrTCTIONS PLOTTED ON THE GRAPH» MAXIMUM GF 10 
MAX = NiUMBER GE DATA PDINITS IN EACH FUNCTIO'N 
BETA = SCALE GiF THE Y AXIS (ABSCISSA) IN UNITS OF TIME 
BU AND @L = THE UPPER AND LOWER LIMITS OF THE X AXIS (OiRDINATE). 

IF BUi=BL=: 0 ,a THE THE MAXiHUM AND MINIMUM VALUES OF 
G THE FUINCTION ARE USED AS THE LIMITS, 

e NZERG = 1 WILL POT A REFERENCE LINE AT X=0. IF A REFERENCE 

e LINE IS NOT DESIRED SEt NZERG = 0. 

e 

INTEGER BLANK, DASH, DOT* DOLAR, AM(10> 

D I ME N S I G N K < ( 10 ) , LINE ( 101) 

DIMENSION FNCT(5,500) ,TIME<BaO) 

DATA NiPR / 6 / 

DATA BL ANK , D A SH ,D OT , A M ( 1 ) , A M ( 2 ) , A Ht 3 ) , AM « A ) , AM (5 ) , AM ( £> ) , A M ( 7 ) « PLOT* 

1 AM( 8) ,AM(9) ,AM( 10) / • •, •-•, •!», *2*, *3*, M*, 'B*, PLOT* 

2 *6*, *7*, *6*, »9», «0* / 

DATA DOLAR / 'S.* / 

IF (BU .EQ. 0.0 -AMD. BL .EG. 0.0) GG TO 210 
BIG = BU 
SMALL = DC 
00 TO 507 
210 8IG=FMCT(l,l) 

SMALL=BIG 

DO 212 M = 1, W:p 

DO 212 L = 2, MAX 
IF (FNiCT(M,L)-BIG> 260,260,213 
213 BIG=FNeT(MtC) 

260 IF (FNiCTCM,L)-SMALL) 214,212,212 
21A SMALL = FfJCT(M,L) 

212 eONTINiUE 
5 07 peak=big-smal:l 

apeak = ABS(PEAK) 

AMlD=(BIG+SMALL)/2.a 


seALE=APEAK/80.0 

UiRITE CMPR»112> AHID» BIB* SMALiLt SCALE 

112 FORMAT H Ml tAX, *l -CUR VE 0 R© I^A TE • * 5X , • HID«VALUE = ^ ♦ 

2 E 1 0 . 3 • 2 1 X f * M A X - V A L U E = • *E 1 0: . 3 1 8 X * • M I N -V A L UE = • t 

3 E 1 0 . 3 . / / * 2 7 X » • S C A L E = * t E 1 0 . 3 , 2 X , * P ER L I WE * t / / ) 

WRITE (N'PR,li3) BETA 

113 FOiRWAT UH , AX i • AB SC ISS A * , 13X r * SC A LE = » » 

2£ 10,3»2X, * PER L I N E * » 62 X , • W -R P S * » // ) 

IF C APEAK .NiE. 0.0 ) GOTO 3 5 
WR ITE {NFR t30 ) 

30 FORMAT (lH0»AXt******* ******************** ******<^****************** 
1.********»****/15X, *Ni0 PLOT WILL BE PRODUCED HiERE BECAUSE*/15X» 
2»FUNCTI0iN MAX VALUE = FUiNCTlOW MIN VALUE. */15Xf * THE FUNCTION IS CO 
3NSTANT. TO RECEIVE A *715X,»P10T OF THIS FUINCTION, PREDETERMINED* 
A/15Xr*UPPER AN® LOWER LIMITS BU AN© BiL MiUST*/15X,*eE SPEC IFI ED . • /5 

6***) 

GOTO A I 
35 CONTINUE 

DO 38 I = 1,101 
39 LIN1E( I)=®ASM 
LINE(li) r DOT 
LINEC51) ? DOT 
LINiElOl) = DOT 
DO A1 d = 1*«AX 
DO AO K=1,NP 

K K ( K ) = 80. 0 * ( F NiC T ( K » d ) - AM 10 ) /APE A K + 51 .0 
KKK = KK(R) 

IF(KKK) 51,51,52 
51 LTNC(l) = QOLAR 
60 TO 5A 

52 IF(KKK-lO'i) 50,50,53 
53 LINE (101) = DOEAR 
GO TO 5A 

50 LINE (KKK) = AM(K) 


ij u U u 


i 


1 


'GOTO 12 

11 ge:hter=o..o 

L€TTER=-1 

12 eO'WTIN'UE 

00 4 KK=1*J 

kkk=kk+j 

YY(KK>=FUSieT<K<»l) 

YY<KKKJ =FUiM.eT(KKf 1 > 

DO 5 L=lfJd 
Y(L)=ftJNGT(KK*L) 

IF m(KK).GT.Y(L))YYCKK)=Y(L) 

IF < YY( KKK) .LT.YCL) ) YY(KKK)=Y(L) 
5 CONTIN'U'E 
4 CONTIWiUE 
d d ' J — 2 *■ d 

CALL SCALE ( YY *SH *ddd »K ) 

CALL SCALE (TIM€» SW* dJ, K) 


i 


X-SCALE OVERRIDE CAW BiE DONE MERE 3Y SPECIFYING TIMECdd + D AS 
TMiE VALUE OF TMIE FIRST TIC MARK ON THE X-AXIS AND TlKE<dd+2) AS 
TMIE CHANGE IN X VALUE BE TWEEN TWO TIC NARKS. 


CALL AXIS (0 .0 »G ENTER f ' • tLETTER *SW • 0*TI«E (dd + 1 ) ,TIN'E *dd + 2 ) ) 
CALL AXIS (0.0,0.0,« • , 0 ♦ SM , 90 . 0 , Y Y ( dd d+1 ) , Y Y t ddd+2 > ) 

DO S NP=l*d 
DO 9 LL=l,dPTS 
Y{LL)=FUiNcf (NIP* LL) 

9 CQNTINUE 

Yldd + l) = YY(ddd“H) 

Y( dd+2)=YY(ddd+2» 

2 CALL LINE (TIME* V, dPTS, 1* 0* 0) 

S eONTINUIE 

CALL PLOT (SU+2.0* 0.0* -3) 

RETURN 
E N,D 



tn 


n o D n o 


S'* CONTINUE 
4 0 eCNTINiUf 

IF { NZERO .fvlE. 1 ) GOTO 44 

KKK = e0.0 * (“AMID)/APEAK ♦ 51.0 

LINE(KKK) = AHUO) 

44 CONTI NiU E 

WRITE (NPRt 39) J* LINEi TIM€<J) 

39 FORMAT HH r4 X , 1 3 « 1 X » 1 0 1 A 1 1 1 X * E 1 0 . 3 ) 

DO 43 L = ItlOl 
43 LINEtL) = BLANK 
L I NE ( 1 ) =D OT 
41 CONTINUE 

WRITE (N.PRt 442 ) 

442 FORMAT ( IH « 26X,***** » •PLOTTING ROUTINE COMPLETED** 12X* •***PLQT* 
1 * ) 

RETURN 

END 


S U 3 R 0 U T I N E M C PLOT < F U NC T * T T M E * d J , J , NO M AG 1 ♦ S W * D A T A ) 

THIS SUB'ROiUTINiE PLOTS GRAPHS ON CAL COMP PLOTTER 
VARIABLES ARE DEFINED AS FOLLOWS: 

FU'NiGT= ARRAY OF Y VALUES 
TIM!E= array of X VALUES 
d= MUMBiE'' OF PLOTS PER GRAPH 
dd= NiUHBER OF POINTS' PER PLOT 
C SU= X AXIS LENGTH IN INCHES 

C SH= Y AXIS LENGTH IN INCHES 

OIMENSICN FUNGT (5,500), TIMIE (500 ), YfSOO), DATA (1024), YY(500) 

SW=10.1 

diPTS = d>d 

SH=7 .5 

K=1 

IF (NOiMAGl) 11, 11, 10 
10 CENTER = SiH 
LETTER-0 


SUBROUTINE G A LC { R < V » OT ♦ X , ID I f REtRDOT iK l»K2f RDT,RVI » L ^ R(« , Y , I DR I ) 

C FILTER 1 BLOCK DIAGRAM 1 
I Nip L I G I T REAL * 0 < D t K * R * V ) 

DIMENSION X(5O,5),Y<25O0) 

IDRiST = IDRI - 50 /ID I 
IDIST = I0I-5O+<5O/IDI)*IDI 
A = GI.Q 

DO 5 0 LAiP=l»IDI 

C ACCELERATION CAN BE ADDCD FOR A PERIOD OF 10*IDI*DT MiERE» L = A1 TG 50 
IFt (L.6T .40 ) .AND.<L.LE.50) ) A=0.0 
V=V+A*DT 

C NEW RANGE » R (TRUE RANGE) 

R=R+V*DT 

e ADD RANDiQH NO.» RM (NOISY RANGE) 

RM^ = R + Y(LAPr' 

G NOISY RANGE MINUS PREVIOUS ESTIMATED RANGE 
RV-RM -RE 

C ESTIMATED RANGE RATE* RD0T 
RQ0T=RD0T+K2 *RV*DT 
G ESTIMATED RANGE, RE 
R0I=RDI+RDOT*DT 
RVI = Rtf I+iKl*RV*DT 
RE=:RDI+RVI 

IF(L .LT. IDRIST) GOTO 50 

IF(L .EQ. IDRIST .AND. LAP .,LE. IDIST) GOTO 50 
LAPX = (L-1)*IDI ♦ LAP - ((I0RTST-1)*IDI ♦ IDIST) 

XtLAPX*!) = RM 
X(LAPX,2) = RE 
X(LAPX,3) = ROOT 
X(LAPX,4) = RE-R 
XtLAPX, 5) = RD'0T-V 
50 CONTINUE 
RETURN 
END 


- » 
i 


-4 


-3 


V 


SU BR oy T I N E e A L G ( R , V , DT t X » I D I » RE » RD OT * K 1 » K2 » R D I » R V I * L ♦ R M * y » I DR n 
e ALPHA -8 ETA FILTER 
e BETA = ALPHA **2/ 12- ALPHA) 

I M P L I C I T R E A L *e C D f K t R » V ) 

D I ME^S I ON X < 50 t5 ) « Y 12 &0 0 ) 
lOiRIST = IDRI - 50/ ID I 
IDIST = IDI-5O+(5O/I0D* lOI 
A =0 . 0 

D O 5 0) LAP = 1,IDI 

C ACCELERATION CAN BE ADDED FOR A PERIOD OF 10.*IDI*DT HERE, L = 41 TO 50 
I F ( ( L . G T . A 0 ) . A R1 D • < L . L E • 5 0 ) ) A = 0 • 0 
V::V + A*DT 

e MEW RANGE, R (TRUE RANGE) 

R=R+V*DT 

C ADD RAND'GM NO., RM (MiOISY RANGE) 

RM = R * Y(LAP) 

RE = RDI*K1*( RH-RDI ) 

RDOT = R00T*K2/DT*(RM-RDI ) 

RDI = RE+DT*RDOT 

IFCL .LT. I DR 1ST) GOTO 5 0 

IF(L .EG. IDRIST .AN®. LAP .LE. IDIST) GOTO 50 
LAPX = (:L-i)*I0I ♦ LAP - ((IORIST-l)*IDI + IDIST) 

XfLAPX,!) = PM 
X(LAPX,2) = RE 
XILAPX,.-?) = ROOT 
X (LAPX, A ) = RE-^R 
X(LAPX,5) = RDOT-V 
50 CONTINUE 
RETURN 
END 


CM 

ro 



SUBROUtINC CA.LG <P»V,DT,X,IDI ,RE,R00TtKl.»K2tRDI»RVl*L*RM,Y»IDPI) 

G FILTER 3 S MG OTHER 

I MP L I C IT R E A L * 8 ( D * K » R » V > 

DIMCNSIGW X<50,5) ,REG (18) »Y(2500) 

IF (L .6T. 1) GOTO 10 
DATA RE 6/18*0.0 00/ 

10 GOWTINiyC 

IDiRIST = IDRI - 50/IDI 
lOIST = IDI-5O*(5O/T0U*IDI 
A = 0 .0 

§0 50 LAiP = l*IDI 

C ACCELERATION GAN BE ADDED FOR A PERIOD OF 10*IDI*DT HERE* L=A1 TO 50 
rFUL.GT.AO ) .ANiD. ( L.LE.50 } )A-G.D 
V=V+A*DT 

C NiEW RANGiE* R (TRUE RANGE) 

R=R+V*DT 

C ADD RANDOM NO.* RM (NOISY RANGE) 

RN = R •* Y(LAP) 

REG(18) = REG(17) 

REG(17) = REG (16) 

REG (16) = REG (15) 

REG (15) = REG(IA) 

REG(IA) = REG (13) 

REG (13) = REG (12) 

REG (12) = REGdl) 

REG(ll) = REG (10) 

REG (10) - RE6( 9) 

R:EG( 9) r RfG( a) 

REG( 8) = REG( 7) 

REG( 7) = REG( 6) 

REG( 6) = REG( 5) 

REGt 5) = REG( A) 

REG( A) = REG( 3) 

REG( 3) = REG( 2) 

REG( 2) = REG( 1) 


cr> 

c»> 




REG( i> = RM 
DS.UM = 0.0 
DO 3 0 I =1,6 

DS.UM = DSDM ♦ REG(I) - REG(l4l2) 

30 CONTINUE 

G ESTIMATED RANGE RATEt RIOT 

RDGTT = DSUIM * A. 5/ (18.0 * OT) 

ROOT = ROOTT/2 3.0 
OS DM. = 0.0 
DO AD 

DSUM = DISUM ♦ REG(I) 

AO CONTINUE 
C ESTIMATED RANGE* RE 

RE = DS'UM/16.0 ♦ ROOT *DT *7.5 
IF(L .EO. 1) RE = RM 
IFtL .EO. 1) root = V 
IF<L .LT. I DR 1ST) GOTO 50 

IF((L .EQ. IDRIST) .AND. (LAP .LE. IDIST)) GOTO 50 
LAPX = (L-1)*IDI ♦ LAP - (( IDiRIST-1 )*IDI + IDlST) 
X(LAPX,1) = RM 
X(LAPX,2> = RE 
XI LAPX *3) = ROOT 
X(LA?X*A) = RE-R 
X(LAPX,5) = RDOT-V 
50 CONTINUE 
RETURN 
END 


4 :. 



SUER OU T I NE C A L C < R » V 1 D T » X , I D I , R E » R 0 0 T ♦ K ] , K2 1 R D I f R V I « L r R M , Y » I D R I ) 
e FILTER A DESIGN 1 

IMP'LIGIT REAL*8(0tKrR,V) 

0J ME MSI ON X ( 5 0 1 5 ) » R E G t 4 ) * Y (2 50 0 ) 

IF (L .6T. 1) GOTO 10 
DATA REG/4+0.GD0/ 

10 GONTlNiUE 

I DR 1ST = I DR I - 5 0/101 

IDIST = IDI-50+t5O/XDI>*IDI ^ 

A =0.0 

DO 50 LAP=1,IDI 

e ACGELERATION GAN BE ADDED FG'R A PERIOD OF 10*IDI*DT HERE, L = 41 TO 50 
IFC (L.GT.40) .AN0.(L.LE.56) ) A=0*0 
V=V+A*DT 

C NEW RANGE, R (TRUE RANGE) 
r=R+V*DT 
RM = R ♦ Y(LAP> 

REG< 4) = RE Of 3) 

REG( 3) = REGC 2) 

REG( 2) = REG( 1) 

;REG( 1) = RM-RE 

DELTA = (REG(l)+REG(2)*RE6t3)*REGC4)>/4.0 

RE = DELTA + RE 

ROOT = DELTA/DT 

IF(L .LT. IDRIST) GOTO 50 

IFtL .E(S* lORIST .AND. LAP ..LE. IDIST) GOTO 50 
LAPX = (L-1)*IDI ♦ LAP - UIDRI ST-1) *101 ♦ IDIST) 

X(LAPX,1) = RN 
X(LAPX,2) = RE 
X (LAPX, 3) = ROOT 
X (LAPX, 4) = RE-R 
X(LAPX,5) = RDOT-V 
5B CONTINUE 
RETURN 

END ^ 

cr» 

cn 



S U© R 0^u T I N E e A L C ( R ♦ V , 0 T , X , I D I » R E » R D G T » K 1 * K 2 » R D I f R V 1 1 L f R K * V f 1 0 R I ) 

C FILTER 5 

C FOiR THIS FILTER, GaLCULATIC-NS ARE MADE EVERY 16 SAMPLES. 

C THIS MEANS THAT FOR EVERY RUN THROUGH GALC IDI/IG VALUES ARE 
e OETiERMINIED RATHER THAN IDI VALUES. TO ALLOW FOR THIS, IDI IS CHANGED 
C TO IDI/16. THIS IS NECESSARY REeAUSE ALL THE VALUES CALCULATFD 
C IN THE LAST RUM THROUGH CALC ARE RETURNED TO BE USED IN VARIANCE 
e eALGULATlONS , 

C 

G *NOTE* - THE IMPLICATION OF THIS CHANGE IS THAT IDI MUST BE AN 
G INTEGRAL MULTIPLE OF 16. ALSO, IDR.I*!©!/ 16 MUST BE GREATER THAN 50. 
e THESE LIMITATIONS APPLY FOR FILTER 5 ONLY. 

C 

IMPilGIT REAL*8 (D,K,R,V) 

DIMENSION X (50,5 ),REG< 16 ),Y( 2500) 

IF (L .GT. 1) GOTO 10 
DATA RE 6/16*0.000/ 

10 GGNTIN'UiE 

IDRIST = IDRI - (50/(101/16)) 

I0I ST = IDI - ( (5 0 - ( 50 / ( I 0 1 /1 6 ) ) * ( I DI/16 ) ) *16 ) 

A = 0 .0 

DO 50 LAP=i,IDI 

C ACCELERATION CAN BE ADiDED FOR A PERIOD OF 10*IDI*DT HERE, L = A1 TO 50 
IF( (L.GT.AO) .ANQ.{L.LE.50))A=0.0 
V=V*A*UT 

e NEW RANGE, R (TRUE RANGE) 

R=R+V*DT 

C ADD RANSOfM NO., RM (NOISY RANGE) 

RM = R + Y(LAP) 

REG (16 ) = REG (15) 

REG(15) = RE6UA) 

REG(IA) = RE6(13) 

REGtl3) = REG(12) 

BEG(12) = REG(ll) 

REG(ll) = REG (ID) 


O 2! 


O 

c 

> 


u 

Q 


a\ 

a\ 




'‘I 






REG (10) 

r 

REGt 

9) 

REG( 

9) 

= 

REG( 

R) 

REGC 

R) 


REG( 

7) 

RE G ( 

7> 


REG( 

6) 

RE G ( 

6) 


REG( 

5) 

R EG ( 

5> 


PEG ( 

4) 

RE G ( 

4 ) 


REG( 

3) 

REG( 

3) 


REG( 

2) 

RlEGt 

2) 

- 

RE6( 

1) 

REG ( 

1 } 


RM-RDI 


IF((LAP“(LAP/1;6>*16) .NC. 0) GOTO 5 0 

3 sun =0.0 

DO 20 I=1,S6 

QSU'M = DSU'M + REGtl) 

20 GQ'NT.IMJiE 

ROOT = DSUM/16.0/DT/i6.0 

R0I = RDI ♦ DSU'M/IG.O 

RE “ RDI + 0 SUM /1 6. 0*7. 5/ 16.0 

IF(L .EQ. 1) RE = RM 

1F<L, .EQ. 1) ROOT = V 

I F(L .LT. TOR I ST) GOTO 5 0 

IFtL .EQ. IDRIST .AND. LAP .LE. IDiSt) GOTO 5(1 
LAPX - <(L-1)*IDI * LAP - ({IDiRIST-l)*IDI + IDIST))/16 


X ( LAPXfl) 


RM 

X (LAPX, 3) 

t 

RE 

X (LAPX ,3) 

z 

RDOT 

X (LAPX ,4 ) 

r 

RE-R 

X (LAPX, 5 ) 

= 

root-v 

continue 



RETURN 



ENiD 




APPENDIX B 


The following pages contain an assembly language listing of 


program FILTER 3. 


FILTER 3 


DECEMBER 1 


♦ 

^ filter 3 


00 1*10 


□PC 

BO 000 

0000 

XTEHP 

priB 

2 

0002 

S©N 

PME 

1 

0 003 

CUT 

PM® 

1 

0 004 

MLTCEM 

PHB 

1 

0 005 

MLTCD2 

PME 

1 

00 OS 

MLTPRl 

PMB 

1 

0007 

NLTPRE 

PMB 

1 

0003 

MLTPR3 

PME 

1 

0009 

PRDIU 

PMB 

1 

OOOfi 

PPDBE 

RMB 

1 

00 OB 

PPQB3 

RMB 

1 

OOOC 

PRafi4 

RMB 

1 

0 0 Oil 

ppni(5 

RMB 

1 

OO0E 

PI • 

RMB 

1 

00 OF 

P2 

RMB 

1 

0010 

P3 

RMB 

1 

00 1 1 

B1 

RMB 

1 

0012 

B2 

RMB 

1 

0013 

B3 

RMB 

1 

0014 

SUMl 

RMB 

1 

0015 

SUM2 

RMB 

1 

00 1 S 
OOIE 

SUMS 

RMB 

□Ri| 

.1 

SOOlg 

OOli 

RTEMP 

RMB 

2 

0020 

0044 

REbsi- 

RM® 

ORi 

3 13 

30044 

0044 

PCNT 

RM® 

1 

0045 

SKPPL© 

RMB 

1 

004S 

PCPCTl 

RMB 

1 

0047 

PCPCTE 

Rf'lB 

1 

0043 

PCPCT3 

RMB 

E+ M p 

1 

r 

U 1 J “ 

OOdfi 

PBOTl 

RMB 

X 

1 

004® 

P00T2 

RMB 

1 

0043 

Ptl 

RMB 

1 

004B 

Pit 

RMB 

1 

004E 

TSUMl 

RMB 

1 

004F 

tSUMt 

RMB 

1 

0 05 0 

T3Ufl3 

RM.? 

1 

0051 

SPMP 

RMB 


0053 

RPiMPfR 

RMB 

*3i 

0 055 

BPTPTR 

RM6 
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FILTER 3 


OeCEMBER 19,. 1973 


0 057 




REPTR 

RMB 

C. 



0 059 




RBDTPT 

RNB 

2 



005B 




CflMPriB 

RMB 

2 



005B 




DIGIT 

RHE 

1 



005E 




TYPE 

RMB 

1 



005F 




ROUND 

RMB 

1 



006 0 




NFLflG 

RMB 

1 



0 0^. 1 




Cl 

RMB 

1 



0062 




C2 

RMB 

1 



0 063 




D1 

RMB 

1 



0064 




B2 

RMB 

1 



0065 




DIFFl 

RMB 

1 



0066 




DIFF2 

RMB 

1 







♦- 









♦ 





0200 





DRG 

$0200 



0200 

C6 

01 


FLT 

LDFl B: 

«f.’I:01 



0202 

Ii7 

5F 



STR B 

RDLINB 



0204 

CE 

00 

01 


LBX 

«.‘I:0G01 



0207 

DF 

51 



STK 

saMP 



0209 

GE 

30 

00 


LBK 

"33 00 0 



02 OG 

DF 

f 



tVA 

REPTR 



02 OE 

GE 

c*y 

00 


LltK 

^J$3@00 



0211 

BF 

59 



STA 

RBQTPT 



0213 

GE 

04 

:3C 


LDA 

«INIT1 



0216 

EB 

El 

4B 


JSR 

:i:E14E 

MIN I BUG 

PBaTRl 

0219 

BB 

El 

33 


JSR 

3E133 

MINIBUG 

INCH 

02 1C 

97 

5E 



:STR a 

TYPE 



02 IE 

44 




L3R a 




02 IF 

24 

05 



BCC 

INPUTl 

' 


0221 

CE 

20 

00 


lb:^ 

jj:I> 2 y 0 U 



0224 

20 

OB 



BRfi 

INPUTS 



0226 

44 



INPUTl 

LSR a 




0227 

24 

05 



BCC 

INPUT2 



0229 

CE 

jE**9,* 

00 


LDK 

"32 300 



022C 

2:0 

03 



BPa 

INPUTS 



022 E 

CE 

12 

00 

INPUT2 

LDK 

"31300 



0231 

BF 

crcr 

•-!' 


INPUT:? 

ft 

DfiTPTR 



0233 

CE 

04 

9E 


LBK 

FINITE 



02:36 

Ml 

El 

4B 


J:3R 

:£EI4B 

MlNIBUG 

pDaTai 

0239 

EB 

EO 

F@ 


J3R 

3E0F3 

MINI BUG 

EaBDR 

023G 

BF 

SB 



;STK 

;S:aNPND 



023'E 

CE 

OG 

£0 


LBK 

"REGS 



0241 

6F 

00 


FLTN 

CLR 

O.-K 



0245 

0© 




INK 




0244 

3C 

00 

44 


CPK 

t^RCYT 



0247 

26 

F8 



EME 

FLTR 



0249 

BE 

55 


FLTl 

LBK 

IiaTPtR 



024B 

EE 

00 



LBK 

U f A 





1 
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1 

FILTER 

3 





DECEMBER 19y 197S 


G£‘‘1D 

BD 

04 

£C 


JSR: 

SHFTR ■ 

1 

0E50 

7F 

0 0 

14 


CLR 

SUMl 

1 

Oc-53 

7F 

00 

15 


CLR 

SUMS 


0£56 

7F 

00 

16 


CLR; 

SUMS 

1 

0S59 

CE 

00 

20 


LEK 

JtREHS 

t 

G£5C: 

C3 

06 



L.DR © 

a S Of 6 


' G25E 

ri7 

44 



STR © 

RCMT 


0£60 

7F 

00 

45 


CLR: 

SKPFL S 

1 

G£t.3 

DF 

0 0 


FLT£ 

STK 

XTEMP 


0£€'5 

CE 

00 

OE 


LE7; 

ijR 1 


. . 0£6S 

fi6 

06 


FLT3 

LDR R 

6 !• r; 


0£fe.fi 

R7 

00 



STR R 

0..X 


0£|“.C 

OS 




IRX 



0£6JIi 

SC 

00 

11 


C-PK 

t^©l 


027 0 

26 

F6 



©RE 

FLT3 

•1 ' 

0£?£ 

HE 

00 



LEK 

KTEMP 

r» 

0£?4 

SC 

00 

£C 


CPH 

{sS00£C 


0£?7 

26 

03 



©ME 

FLT4 

• i. 

- tj27'5» 

7R 

00 

45 


DEC 

SKRFLS 


y w 1** L 

SC 

00 

3S 

FL T4 

CPX 

“S 0 OSS 


0£7F 

£6 

03 



©ME 

FLT5 

55 • 

0£81 

•• f 
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61 


LOdP£ 

STR R 

Cl 

0577 

17 

6£ 



SIR B 

C£ 

0579 

£0 

06 


■ 

ERR 

L00P4 

057E 

EE 

00 


LOOPS 

LIK 

0 j K 

05711 

EE 

00 


LOOPS 

LIK 

OjK 

057F 

IF 

61 



STK 

Cl 

05:;!l 

Bl 

05 

B1 

LOOP 4 

JSR 

CDNV 

05ili4 

IE 

00 



LIK 

KTEMP 

0586 

08 




INK . 


0587 

0© 




IMK 


058® 

IF 

00 



STK 

KTEMP 

05© R 

®C 

00 

5B 


CPK 

isSRMPNl 

058D 

£6 

Cl 



BME 

LODPl 

05:®F 

IE 

53 



LIK 

RRMPTR 

0591 

03 




IMK 


0593 

0® 




INK 


0593 

IF 




STK 

RRMPTR 

0595 

IE 

55 



LIK 

IRTPTP 

0597 

03 




INK 


0598 

0S 




INK 


0599 

IF 

55 



STK 

IRTPTR 

059B 

IE 

57’ 



LIK 

REPTR 

©591 

03 




INK 


05 9 E 

03 




INK 


059F 

IF 

ST' 



STK 

REPTR 

05R 1 

EE 

59 



LEK 

RlOTPT 

05R3 

08 




INK 


05R4 

03 




INK 


05R5 

IF 

59 



STK 

RIOTPT 

05R7 

IE 

51 



LIK 

SRMP 

05R9 

9C 

5B 



Cf-\K 

Sl^MPNl 


MlfiliUG 


OUTCH 


MIHIBUe 


PiiftTRl 


PDRTRl 
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FILTIR 3 





0ECEMEER 

19> 19?8 

05F<B 


05 



0EQ 

PEM0 


05FID 

08 




I NX 



i:i5flE 

0F 

51 



STX 

SFIMP 


05B0 

£0 

91 



0RP 

LOOP 

, 

05££' 

©6 

OF* 


FEN0 

L0fl Ft 

<?©OFl 


05E‘*r 

C6 

05 



LDFi B 

«3:0.5 


0506 

00 

El 

26 

PEN 01 

JSR 

©El £6 

MINI BUG DUTCH 


5R 




0EC 0 



t:i5EFi 

£6 

Ffi 



ENE 

PEMDl 


OSEC 

3F 



♦ 

swr 











0-5Brt 

?F 

00 

60 

CDNV 

CLR 

MFLfiG 


05C 0 

?0 

00 

61 


TCT 

cr 

- 

mscs 

£F 

03 



0PL 

CDNVl 


05C5 

?F 

00 

60 


0EC 

NFLRG 


05CS 

CE 

*?- 
C ( 

1 0 

CQNVl 

L0K 

:f©£?10 

1 0!. 0 00 

OSCE 

00 

05 

E8 


J3R 

Ft=iCTaR 


05CI 

CE 

03 

E© 


L0X 

'»S 03iEB 

1 > 000 

0501 

00 

05 

E8 


J3R 

FfiCTDR 


0504 

CE 

00 

64 


LUX 

if© 0 ij64 

1 00 

050? 

00 

05 

ES 


JS'R 

FRCTOR 


O50t^ 

CE 

00 

OFi 


L0X 

WSOOOFI 

10 

0500 

00 

05 

ES 


J£R 

FFiGTaR 


05E0 

•H 

6£ 



L0FI Fit 

CS 


0512 

©FI 

30 



ORF) FI 

«:f30 


0514 

00 

El 

£6 


JSR 

SEle6 

MIN I BUG DUTCH 

051? 

©S' 



♦ 

RT3 







♦ 


• 


05E© 

BF 

63 


FACTEIR 

STX 

01 


05If^ 

?F 

00 

50 


CLR 

BISIT 

‘ 

■ 050 

00 

06 

00 

FFICT 

JSR 

SUBe'Ey 


05F0 

£0 

OC 



EMI 

FRCTN0 


05Fc; 

7F 

00 

60 


CLR 

NFLFIG 


05F5 

7C 

00 

50 

FFICTI 

INC 

01 GIT 


05F© 

0E 

65 



L0X 

01 FFI 


05FFI 

0F 

6 1 



STX 

Cl 


05FC 

£0 

EF 



ERFI 

FRCT 


05FI 

?B 

00 

6 Cl 

FF'CTMB 

1ST 

NFLFIS 


06 01 

ES 

F2 



EMI 

FFCTI 


06 03 

©6 

51) 



L 0Ff Ft 

01 SIT 


0605 

©Ft 

30 



DRF) Ft 

is$30 


060? 

SB 

El 

26 


JSR 

©El £6 


06GFi 

39 



•# 

RTS 




♦ S rsVTt SYliTRfiCTiari ROUT INI 


SUlf:£3V L0?2| B Cg 


iIigOB 06 62 
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FILTER 3 


DECEMBER 19.» 197S 


O'iori 

96 

61 


LDR 

R 

Cl 

06 OF 

DO 

64 


SUE 

E 

D2 

0611 

92 

63 



R 

III 

0613 

D7 

66 


STR 

E 

DIFF2 

0615 

97 

65 


SIR 

R 

DIFFi 

0617 

39 



RTS 






♦ 







* 




0613 

OD 


HERD! 

FCE 


3 0 D 3 OR • 3 OR 5 3 OR j 3 OR 

0619 

OR 

OR 





06 IB 

OR 

OR 





06 1 Ii 

20 



FCC 


.••• SRMPLE ! RRMP 

U6 IE 

53 

41 





0630 

4D 

50 





0632 

4C 

45 





0624 

20 

20 





0626 

20 

20 





0623 

21 

20 





062R 

20 

20 





OSSL’ 

20 

52 





0621 

41 

4B 




. 

0630 

50 

20 





0632 

20 

20 





0634 

20 

21 





0636 

20 






0637 

04 



FCE 


304 

0633 

20 


RHERD 

FCC 


^ RRMP 

0639 

20 

20 





U63B 

52 

41 





063B 

4D 

50 





iji62*F 

20 

£0 





0641 

20 






0642 

04 



FCE 


304 

0643 

20 


HHERK 

FCC 


NOISE / 

0644 

20 

20 





0646 

41 

4F 





0643 

49 

53 

• 




0641^ 

45 

20 





064 C 

20 






06 4B 

04 



FCE 


304 

064E 

52 


RMHERD 

-FOC 


.-•■RRNP+NOISE-' 

06 4F 

41 

40 





0651 

50 

23 





0653 

4E 

4F 





0655 

49 

53 





0657 

45 






0658 

04 



FCE 


304 

06S9 

ao 


N:ERf(2 

FCC 


! RE ! RDOT^- 

06S8 

fi 

go 
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FILTER 3 BECEMBER 19, 1973 


G35C 

20 

SO 




065E 

c'O 

20 




0660 

52 

45 




0662 

20 

U 




0364 

20 

SO 




0666 

SO 

21 




0363 

go 

so 




066fi 

go 

20 




066C 

5g 

44 




066E 

4F 

54 




0670 

OB 


HERB 3 

FCE 

SOB, 3 OR, 304 

0671 

OR 

04 




0673 

2B 


HERBS 

FCC 


0674 

2B 


HERD4 

FCC 



0675 

2B 

2B 




0677 

2D 

2B 




0679 

gn 

2D 




06 7B 

2B 

2B 




067B 

gB 

£B 




067F 

04 



FCB 

*04 

0630 

gl 


HERDS 

FCE 

*21**04 

0631 

04 





0632 

00 

OB 

HERD6 

FBE 

*0D,*0R,S20, 

0634 

00 

OR 




0636 

00 

20 




0633 

00 

20 




U 63 F 1 

00 

04 




0633 

go 


HERB7 

FCC 

/ 1 / 

1 II 6 SD 

20 

£6 




063F 

20 

21 



' 

0691 

20 

20 




0693 

20 





0694 

04 



FCB 

*04 

0695 

©3 


MINUS 

FCB 

*03, *21, *04 

0696 

2B 

04 








END 



HQ 

1 ERRQR 

c;s> DEI 

rECTED 



SVMBQ 

L 'RBLi: 





R1 

00OE 

R2 

ij'U yF 

R3 

001 0 

RDIi 

0406 

B1 

© 0 1 1 

IS 

0012 

B3 

0013 

Cl 

006 1 

C2 

0062 

CTIT 

OO-03 

CDNV 

05 BB 

COHVl 

05C3 

B1 

0063 

B2 

0064 

BRTRTR 

0055. 

BIFF I 

0065 

DI-'Fg 

0066 

BUS IT 

0 05B 

Li-irt 

0361 

FRCT 

05EB 

FRCTl 

osrs 

FRCTNB 

05FE 

FRCTQR 

OSES 

FLT 

0200 



FILTER 


o 


FLTI 

i:i£4S 

FLTl O 

FLT14 

03£F 

FLTS 

FLT4 

Ci£7C: 

FLT5 

FLT? 

0£ft 0 

FLTR 

HEAIfc' 

0®59 

HERrfS 

HERDS 

06 SO 

HERD® 

HERD® 

0670 

INITl 

IMPUTl 

):i££6 

IMPUTE 

L.OUP 

0543 

LDdPl 

LQPP3 

057B 

LDGP4 

MINUS 

0695 

MLTC~il 

MLTPRl 

0006 

MLTPR2 

MULT 

0S64 

MULT 1 

MULT3 

039R 

NULT31 

MULT5 

03D5 

MULTND 

NESg 

0405 

NE33 

NEURTE 

03FE 

NFLfi® 

□UT 

0337 

DUTl 

PENDl 

iH5B$ 

PR INTI 

PRINTS 

05£5 

PRINT4 

PRINTR 

04E4 

PRQDl 

PRODS 

00 OB 

PROD 4 

RRMPTR 

0053 

RCNT 

RCRC.T2 

0047 

RCRCT3 

RDOTl 

0O4R 

RDOTE 

REl 

004C 

REE 

REPTR 

0057 

RMERD 

ROUND 

005F 

ROUND 1 

RTEMP 

OOlE 

SRMP 

SSN 

0003 

SNFTR 

SKPFL© 

0045 

SUB 

SIJMl 

0014 

SLIME 

TSUMl 

0O4E 

TSLIME 

TYPE 

0O5E 

KTt'MP 



DECEMBER 19» 

1973 

0EE9 

FL TIE 

OEEE 

0E63 

FLT3 

0E63 

0E34 

FLT6 

0E9D 

0E41 

HERDl 

0613 

0673 

HERXf4 

0674 

06SE 

HERD? 

063C 

043C 

I MITE 

049E 

OEEE 

INPUTS 

0E31 

0550 

LOOPS 

0573 

0531 

LOOPS 

057D 

0004 

MLTCDE 

0005 

0 00? 

MLTPR3 

00 OS 

03 7C: 

MUL TE 

0394 

03RR 

MULT4 

03B4 

03 F3 

NE31 

0404 

03F9 

NEG4 

03FF 

0060 

NHERD 

0643 

0341 

PEND 

05ES 

05 OE 

PR I NTE 

051 G 

0534 

PRINTS 

053D 

0009 

PROBE 

00 0R 

O0OG 

PRODS 

0 0 0D 

0044 

RCRCTl 

0^046 

0048 

RCRC.T4 

0-049 

IK04B 

RDDTPT 

0059 

004D 

RE3S 

00E0: 

0633 

RNHERii 

064E 

03B!| 

ROUNDE 

03CF 

0051 

SRMPND 

005B 

04EC 

SHFTRl 

0431 

0419 

SUBEEY 

060® 

0015 

SUMS 

0016 

0 04F 
00 00 

TSUM3 

0 050 


APPENDIX C 


The following pages contain an assembly language listing of 


program FILTER 5. 
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FILTER 5 


DECEfIBER 19> 1978 


♦ FILTER 5 


0000 

0000 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 
00 OF) 
00 OB 
0OOC 

00 on 

OO0E 

fiOOF 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
001 E 
001 E 
0020 
0040 
0044 

0044 

0045 
0'046 
0047 
0043 
0049 
0043 
0O4B 
0O4e 
0040 
0LTt4E 
004F 
0051 
0053 
0055 
0057 


XTEMP 

□R6 

RMB 

30000 

2 

SGN 

RMB 

1 

CMT 

RflB 

1 

MLTCDl 

RMB 

1 

MLTcna 

RMB 

1 

MLTPRl 

RMB 

1 

MLTPR2 

RMB 

1 

MLTPR3 

RMB 

1 

PROD! 

RMB 

1 

PRQB2 

RMB 

1 

PR0D3 

RMS 

1 

PR014 

RMB 

1 

pRons 

RMB 

1 

PI 

RMB 

1 

B2 

RMB 

1 

m 

RMB 

1 

El 

RMB 

1 

E£ 

RMB 

1 

B3 

RMB 

1 

sum 

RMB 

1 

SUM2 

RMB - 

1 

SUM3 

RMB 

1 

RTEMP 

QRS 

RMB 

3001 E 
2 

REGS 

RMB 

32 

REGSHB 

RMB 

4 

RCNT 

SR5 

RMB 

30044 

1 

SKPFLG 

RMB 

1 

RCRCTl 

RMB 

1 

RCRCTS 

RMB 

1 

R©QT1 

RMB 

1 

R1QT2 

RMB 

1 

REl 

RMB 

1 

RE2 

RMB 

1 

TSUMl 

RMB 

1 

TSUN2 

RMB 

1 

Tsuri3 

RMB 

1 

SPMP 

RMB 

2 

RPMPTR 

RMB 

2 

DPTPTR 

RMB 

2 

REPTR 

RMB 

2 


RHaTPT RMB 2 


FILTER 5 


DECEMBER 19> 1-??S 


0059 
005B 
00-5C 
005D 
O05E 
005F 

0060 
0061 
006S 

0063 , 

0064 

0065 
0067 


0200 

0200 C6 01 
0202 D7 5D 
0204 CE ‘00 01 
0207 DF 4F 
0209 CE 30 00 
020C DF 55 
020E CE 38 00 
0211 IF 57 
0213 CE 04 IF 
0216 BD El 4B 
0219 BD El 33 
02 1C 97 5C 
02 IE 44 
02 IF 24 05 
0221 CE 20 00 
0224 20 OB 
02c*6 44 
0227 24 05 
0229 CE 28 00 
022C 20 03 
022E CE 18 00 
0231 DF 53 
0233 CE 04 81 
0236 BD El 4B' 
0239 BD E0 F§ 
023C DF 59 
023E CE 0'0 20 
0241 6F 00 

0243 08 

0244 SC 00 44 
0247 26 F8 
0249 7F 00 67 
024C 06 10 


RUB 

2 

RNB 

1 

RME 

1 

RNB 

1 

RMB 

1 

RMB 

1 

RMB 

1 

RMB 

1 

RMB 

1 

RMB 

1 

RME 

1 

RME 

2 

RMB 

1 


□RG 


30200 

LDFt 

B 

,'?301 

STB 

B 

ROUND 

LDX 


«$O'0O1 

STX 


SBMP 

LD^^ 


«33000 

STX 


REPTR 

LDX 


«83800 

STX 


RDDTPT 

LDX 


iilNITl 

JSR 


3E14B 

JSR 


3E133 

STB 

B 

TYPE 

LSR 

B 


BCG 


INPUTl 

LDX 


«32000 

ERB 


INPUTS 

LSR 

B 


BCG 


INPUT2 

LDX 


is:J28O0 

BRB 


INPUTS 

LDX 


«3180O 

STX 


DBTPTR 

LDX 


«sINIT2 

JSR 


3E14B 

JSR 


3E0F© 

STX 


SBMPND 

LDX 

CLR 


UREGS 
Oj X 

I NX 
CPX 


ssRCNT 

BNE 


FLTB 

CLR 


ENiDFLG 

LDB 

B 

is310 


SfiMPMD 

DIGIT 

TYPE 

ROUMD 

NFLPG 

Cl 

C2 

D1 

D2 

DIFFl 

DIFF2 

RDI 

EMDFLG 


FLT 


IMPUTl 


IMPUT2 

INPUT3 


FLTP 


MIN I BUG PDRTRl 
MIN I BUG INCH 


NINIBUG PDHTPl 
MIN I BUG BHDDR 
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FILTER 5 


DECEMBER f?» 1973 


024E 

D7 

44 



STR B 

RCNT 

0£50 

DE 

53 



LDX 

DBTPTR 

0352 

EE 

00 



LBK 

o.*x 

0254 

IF 

65 



STX 

RDI 

0256 

DE 

53 


FLTl 

LDX 

DBTPTR 

0253 

EE 

00 



LJJK 

o.*x 

025R 

DF 

OF 



SIX 

B2 

025C 

DE 

65 



LDK 

RDI 

025E 

DF 

12 



STX 

B2 

0260 

7F 

00 

OE 


CLR 

B1 

0263 

7F 

00 

11 


CLR 

11 

0266 

ID 

03 

FC 


JSR 

SUB 

0269 

DE 

15 



LDX 

SUM2 

026B 

ID 

04 

OF 


JSR 

SMFTR 

026E 

DE 

53 



LDK 

DBTPTR 

0270 

03 




INX 


0271 

03 




I NX 


0272 

DF 

53 



STX 

DBTPTR 

0274 

DE 

’4F 



LDX 

SBMp 

0276 

03 




INX 


0277 

DF 

4F 



SIX 

SBMP 

0279 

9G 

59 



CPX 

SBMPND 

027E 

26 

03 



BNE 

FLTll 

0270 

7C 

00 

67 


INC 

ENDFLS 

0280 

73 

00 

44 

FLTll 

DEC 

RCNT 

0233 

26 

D1 



BNE 

FLTl 

0235 

C6 

10 



LBR B 

»3l 0 

0237 

17 

44 



STB 1 

RCNT 

0239 

7F 

00 

14 


CLR 

sum 

023C 

7F 

0'0 

15 


CLR 

SUM2 

023F 

7F 

00 

16 


CLR 

SUM3 

0292 

CE 

00 

20 


LDX 

«RESS 

0295 

DF 

00 


FLT2 

SIX 

XTEMP 

0297 

CE 

00 

OE 


LDX 

«B1 

0293 

36 

06 


FLT3 

LDR B 

6.0 X 

0290 

37 

00 



STB B 

0> X 

029E 

03 




INX 


029F 

3C 

00 

11 


CPX 

OBI 

0232 


F6 



BNE 

FLT3 

0234 

7F 

O'0 

11' 


CLR 

11 

02fi7 

DE 

00 



LDX 

XTEMP 

0239 

m 

00 



LDX 

0*X 

0gRB 

m 

03 



BPL 

FLT4 

0231 

73 

0'0 

11 


DEC 

11 

0210 

IF 

13 


FLT4 

3TX 

12 

0212 

ID 

03 

E9 


JSR 

BID 

0215 

DE 

00 



LDX 

XTEMP 

0217 

03 




INX 


0213 

03 




INX 
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FILTER 5 


0£E9 

op 

o'-- 

00 

40 

0£EC 

£6 

D7 


0£BE 

96 

14 


0£C0 

97 

4C 


0£C£ 

97 

06 


0£C4 

BE 

15 


02C6 

3F 

4B 


0£CS 

BF 

07 


0£CR 

CE 

06 

40 

0£CD 

DF 

04 


0£CF 

BD 

03 

47 

0£B£ 

BE 

OR 


0£D4 

BF 



0£B6 

B6 

4C 


0£DS 

B7 

06 


0£Dft 

BE 

4B 


02DC 

BF 

07 


OSDE 

CE 

10 

00 

0£E1 

BF 

‘04 


02E3 

BB 

03 

47 

0£E6 

96 

09 


0£ES 

97 

11 


0£EFl 

BE 

OR 


Cf£EC 

BF 

1£ 


0£EE 

BE 

65 


0£F0 

BF 

OF 


0£F£ 

7F 

00 

OE 

iH£F5 

EB 

03 

E9 

0£FS 

BE 

15 


iIigFff 

BF 

65 


02FG 

B6 

4C 


iIigFE 

B7 

06 


03 00 

BE 

4B 


030£ 

BF 

07 


0304 

CE 

17 

£0 

0307 

BF 

04 


03 09 

EB 

03 

47 

030C 

5*6 

09 


0301 

97 

11 


0310 

EE 

OR 


03 1 £ 

BF 

1£ 


03 1 4 

EB 

03 

E9 

0317 

BE 

15 


0319 

BF 

4R 


0311! 

BE 

55 


031 D 

96 

4R 


03 IF 

B6 

4B 


03£1 

Fi7 

00 


03£3 

E7 

01 



BECEMEER 


CPX 


«REGSHD 

BNE 


FLT2 

LDR 

R 

sum 

STR 

R 

TSLIMl 

-STR 

R 

MLTPRl 

LBX 


SUM2 

STX 


TSUM2 

srx 


MLTPR2 

LB^: 


«:f:0640 

STX 


MLTCBl 

J-SR 


MULT 

LBX 


PRI3D2 

STX 


RBOTl 

LBR 

B 

Tsum 

STR 

B 

MLTPRl 

LBX 


TSUM2 

STX 


MLTPR2 

LBX 


«31000 

STX 


MLTCBl 

JSR 


MULT 

LBR 

R 

PRDBl 

STR 

R 

B1 

LBX 


PR0B2 

STX 


E2 

LBX 


RBI 

STX 


R2 

CLR 


R1 

JSR 


RBB- 

LBX 


SUM2 

STX 


RBI 

LBR 

B 

TSUMi 

STR 

B 

MLTPRl 

LBX 


TSUM2 

STX 


MLTPR2 

LBX 


^^31730 

STX 


MLTCBl 

JSR 


MULT 

LBR 

R 

PROIii 

STR 

R 

El 

LBX 


PRDB2 

STX 


E2 

J*iR 


RDB 

LBX 


SUM2 

STX 

I TTiV 


REl 

pep TP 

LBR 

R 

REl 

LBR 

B 

RE£ 

STR 

R 

O.'X 


STFi B 1»X 


r?!. TS*73 


□UT 
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FILTER 5 






0335 

03 




I NX 


0326 

03 




I NX 


0327 

DF 

55 



3 TX 


0323 

DE 

er? 

%-H 


OLITl 

LDX 


032 B 

96 

43 



LDR 

R 

032 D 

B 6 

49 



LDR 

B 

0 ? 2 F 

R 7 

GO 



STR 

R 

0331 

E 7 

01 



STR 

B 

0333 

03 




INX 


0334 

03 




I NX 


iioc'c; 

DF 

57 



STX 


0337 

7 D 

00 

67 


TST 


033 R 

26 

03 



ENE 


033 C 

7 E 

02 

56 


JMP 


033 F 

BE 

4 F 


END 

LDX 


0341 

09 




DEX 


0342 

DF 

59 



STX 


0344 

7 E. 

04 

C 7 


JMP 













♦ 3 

BYTE X 

2 





♦ 



0347 

DF 

GO 


MULT 

STX 


0349 

CE 

00 

02 


LDX 


034 C 

6 F 

00 



GLR 


034 E 

6 C 

00 



INC 


0350 

R 6 

02 



LDR 

R 

0352 

2 R 

OB 



EPL 


0354 

E 6 

03 



LDR 

B 

0 c <36 

BD 

03 

El 


JSR 


0359 

R 7 

02 



STR 

fl 

035 B 

E 7 

03 

. 


StR 

B 

035 B 

6 R 

00 



DEC 


035 F 

R 6 

04 


MULTI 

LDR 

fi 

0361 

2 R 

14 



BPL 


0363 

R 6 

05 



LDR 

R 

0365 

E 6 

06 



LDR 

B 

0367 

BD 

03 

21 


JSR 


036 R 

R 7 

05 



STR 

R 

036 C 

E 7 

06 



STR 

B 

036 E 

R 6 

04 



LDR 

R 

0370 

ID 

03 

12 


JSR 


0373 

R 7 

04 



STR 

R 

0375 

6 R 

GO 



DEG 


0377 

36 i 

E 7 


MULT 2 

LDR 

R 

0379 

R 7 

01 



STR 

R 

037 B 

4 F 




CLR 

R 

037 C 

5 F 




CLR 

B 

037 D 

6 C 

01 


MULT 3 

INC 


037 F 

27 

16 



EEQ 



DEGEMBgR 1 S*.« 19?8 


REPTR 

RDQTPT 

RDQTl 

RDOTS 

OiiX 

1»X 


RDOTPT 

EHIiFLS 

END 

FLTl 

SRMP 

SRMPND 

PRINTR 

BYTE GOMPLEMEMt MULTIPLY 

XTEMP 

0»X 

OfX 

2 >X 

MULTI 

3 jX 

MESRTE 
2 jX 
3 jX 
0 »X 
4s X 
MULTS 
5sX 
6 s X 

MESRTE 
5s X 
6sX 
4s X 
ME34 
4s X 
0 s X 
«:tE7 
IsX 


IsX 

MULT 4 



FILTER 5 


DECEMBER 1 


03:51 

64 

04 



LSR 


4,X 

0333 

66 

05 



RDR 


5,X 

0335 

66 

06 



ROR 


6,X 

03:3? 

£4 

04 



ECC 


MULT31 

033’? 

EE 

03 



BED 

E 

3!»X 

033 B 

B? 

OS 



BDC 

B 

£>X 

033D 

44 



MULT31 

LSR 

B 


033E 

56 




ROR 

B 


03SF 

66 

09 



ROR 


9>X 

03’51 

6-6 

OB 



ROR 


10.-X 

0333 

66 

OB 



ROR 


lliiX 

0335 

SO 

E6 



ERB 


MULT3 

03’5? 

B? 

07 


MULT4 

STB 

B 

7jX 

03 '5 3 

E? 

Oy 



STB 

B 

8yX 

03‘?E 

?D 

00 

5D 

RDUHDl 

TST 


ROUND 

03'3E 

S? 

18 



EEQ 


MULT5 

033.0 

71 

00 

OC 


TST 


PRaD4 

03PI3 

m 

01 



BPL 


ROUNDS 

'J3B5 

7G’ 

00 

OE 


INC 


PR0D3 

03B:5 

S6 

03 



ENE 


ROUNDS 

03BB 

?C 

00 

OB 


INC 


PRODS 

03BH 

£6 

03 



ENE 


ROUNDS 

03BF 

?C 

00 

09 


INC 


PRDDl 

03B3 

?F 

00 

OC 

ROUNDS 

CLR 


PR0D4 

03B5 

7F 

00 

0D 


CLR 


PRQD5 

03B3 

66 

00 


MULT5 

ROR 


0»X 

03BR 

£5 

ID 



ECS 


MULTNl 

03BG 

B6 

OB 



LDB 

fi 

IOjX 

03EE 

E6 

OB 



LDB 

B 

H jX 

03CO 

ED 

03 

El 


JSR 


NE6BTE 

03C3 

B7 

OB 



STB 

B 

10,-X 

03C5 

E7 

OE 



STB 

£ 

11»X 

03C? 

B6 

08 



LDB 

B 

3»X 

03C? 

E6 

09 



LDB 

E 

9).X 

OSCE 

ED 

03 

DC 


JSR 


NES3 

03C1 

B? 

03 



STB 

B 

S»X 

0310 

E? 

0*9 



STB 

B 

9jX 

0313 

B6 

07 



LDB 

B 

?!.X 

0314 

ED 

03 

ES 


JSR 


NE34 

031? 

B7 

07 

- 


STB 

B 

7fH 

0:|M9 

DE 

O'O 


MULTMD 

LDJ< 


XTEMP 

03DE 

39 




RTS 







♦ 








♦ S'S 1 

COMPLEMENT ROUT I NE 





♦ 




03BC 

£4 

03 


NE33 

ECC 


NEeBTE 

0311 

53 




COM 

£ 


OSIF 

SO 

06 



ERB 


NESl 

03E1 

50 



HESBTE 

NE3 

E 
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FILTER 5 


DECEMBER 19.- 1978 


03E2 

ES 

03 

NE64 

ECS 

NEGl 

03E4 

40 



MEG R 


03E5 

20 

01 


ERR 

NEG2 

03E7 

43 


NEGl 

COM R 


03ES 

39 


NE62 

RTS 





♦ 






♦ 3 

BYTE RDDITION ROUTINE 




♦ 



03E9 

D6 

10 

RDD 

LDR E 

R3 

03EB 

98 

OF 


LDR R 

R2 

03ED 

DB 

13 


RDD B 

E3 

03EF 

99 

12 


RDC R 

B2 

03F1 

D7 

16 


STR B 

SUM3 

03F3 

97 

15 


STR H 

SUM2 

03F5 

98 

OE 


LDR R 

R1 

03F7 

99 

11 


RDC R 

B1 

03F9 

97 

14 


STR R 

SUM! 

03FE 

39 



RTS 





♦ 






♦ 3 

BYTE SUBTRACT I ON ROUTINE 







G3FC 

D6 

10 

SUB 

LDR B 

R3 

03FE 

98 

OF 


LDR R 

R2 

0400 

DO 

13 


SUB E 

B3 

0402 

92 

12 


SBC fi 

B2 

0404 

D7 

16 


STR E 

SUMS 

04Cii* 

97 

15 


STR R 

SUM2 

0403 

98 

OE 


LDR R 

R1 . 

04 OR 

92 

11 


SBC R 

El 

04 OC 

97 

14 


STR R 

SUMl 

04 OE 

39 



RTS 



♦ ROUTINE TO SHIFT DHTfi IN REGS 


04 OF 

BF 

IE 


SHFTR 

STX 

0411 

CE 

GO 

41 


LDX 

0414 

E6 

00 


SHFTRl 

LDR B 

041 6 

E7 

02 



STR B 

0413 

09 




DEX 

0419 

3C 

00 

ID 


CPX 

04 1C 

26 

F6 



BNE 

04 IE 

39 




RTS 

04 IF 
0420 

OD 

OR 

OR 


INITl 

FCB 

0422 

0424 

0426 

OR 

OR 

46 

OR 

OR 



Fce 

0427 

49 

4C 





RTEMP 
t;:J 15041 
0>X 
£jX 

«J:001D 

SHFTRl 


$D.* SHii -Eh j EH,* 3'8 


-'FILTER 5 MICROFROeESSOR-'- 
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FILTER 5 DECEMBER 19? 1973 


0429 

54 

45 



042B 

52 

20 



04 2D 

35 

20 



042F 

4D 

49 



0431 

43 

52 



0433 

4F 

50 



0435 

52 

4F 



0437 

43 

45 



0439 

53 

53 



043E 

4F 

52 



043D 

20 


FCC 

SIMULftTIDM RUM/ 

043E 

53 

49 



0440 

4D 

55 



0442 

4C 

41 


' 

0444 

54 

49 



0446 

4F 

4E 



0443 

20 

52 



044ft 

55 

4E 



044C 

0D‘ 


FCB 

3D? 3 ft? 3ft 

044D 

Oft 

Oft 



044F 

49 


FCC 

/irjPUT CHQICE?/ 

0450 

4E 

50 

' 


0452 

55 

54 



0454 

20 

43 



0456 

48 

4F 


■ 

0453 

49 

43 



045fl 

45 

3F 



045C 

OD 


FCB 

3D?.3ft 

045D 

Oft 




045E 

23 


FCC 

/ <SI6MftL=0? M0ISE=1 ? / 

045F 

53 

49 




0461 47 4E 
0463 41 4C 
0465 3D 30 
0467 SC 4E 

0469 4F 49 
046D 53 45 
046D 3D 31 
046F i*€ 

0470 53 • FCC -^SI@HRL+MQISE-3> 

047 1 49 47 
0473 4E 41 
0475 4C 2E 


0477 

41 

4F 




V'T ♦ 



047® 

45 

30 

047D 

32 

29 

047F 

20 


0430 

04 



FCB 


04 


riLTER 5 


DECEMBER 19.* 1973 


0431 

OD 


IMIT2 

FCB 

0432 

OA 

OA 



0434 

4E 



FCC 

0435 

55 

4D 



0437 

42 

45 



0439 

52 

£0 



043B 

4F 

46 



0430 

£0 




04SE 

53 



FCC 

043F 

41 

4D 



0491 

50 

4C 



0493 

45 

53 



0495 

20 

54 



0497 

4F 

20 



0499 

42 

45 



049E 

20 

52 



0490 

55 

4E 



049F 

3F 




04 AO 

OD 



FCB 

04A1 

OA 




04A2 

£8 



FCC 

04A3 

45 

4E 


t 

04A5 

54 

45 



04A7 

52 

£0 



G4A9 

41 

£0 



04AB 

34 

20 



04AD 

44 

49 



04AF 

47 

49 

* 


04B1 

54 

£0 



04B3 

43 



FCC 

04B4 

45 

58 


, 

04B3 

49 

44 



04BS 

45 

43 



04BA 

49 

4D 



04BC 

41 

4C 

• 


04BE 

20 

56 



O4C0 

41 

4C 



04C2 

55 

45 



04C4 

29 

20 



04C6 

04 



FCB 


3D*;!;a.*3i=i 
/NUMBER OF / 


/SAMPLES TO BE RUM?/ 




/<ENTER A 4 DISIT / 


/HEX I DECIMAL VALUE> / 


04 


♦ 


♦ FILTER DATA PRINTER PROGRAM 


♦ 


04C7 

CE 0 0 

10 

PRINTR LDX 

00 1 0 

04CA 

DF 4F 


SIX 

SAMP 

04CC 

CE 13 

IE 

LDX 

it:l:l31E 

04CF 

DF 51 


STX 

RAMPTR 
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FILTER S 






DECEMBER 

13j 1373 

04 D1 

CE 

30 

00 


LDX 


«33000 


04D4 

DF 

55 



STK 


REPTR 


04D6 

CE 

•-*o 

00 


LDX 


00 


04D^ 

DF 

57 



STX 


RDDTPT 


04DB 

CE 

06 

11 


LIii< 


-HEPDl 


04IiE 

£D 

El 

4B 


J$R 


3E14B 

MINIBUS PDPTPl 

04 El 

36 

SC 



LDP 

P 

TYPE 


04E3 

44 




LSR 

P 



04E4 

S4 

OB 



BCC 


PR INTI 


04E6 

CE 

06 

3C 


LD?< 


WNHERD 


04E9 

BD 

El 

4B 


JSR 


:l:E14E 

MINIBUS PDPTPl 

04EC 

CE 

£0 

IE 


LDX 


»3£01E 


04EF 

£0 

17 



ERP 


PRINTS 


04F1 

44 



PRINTl 

LSR 

P 



G4F£ 

£4 

OB 



BCC 


PRINTS 


04F4 

CE 

06 

47 


LDK 


ssRNHEPD 


04F7 

BD 

El 

4B 


JSR 


8E14B 

MINIBUS PDPTPl 

04Ff=l 

CE. 

£8 

IE 


LDK 


«8£81E 


04FD 

SO 

09 



ERR 


PRINT3 


04FF 

CE 

06 

31 

PRiriT£ 

LDK 


^JRHEPD 


05 OS 

BD 

El 

4B 


JSR 


3E14B 

NIMBUS PDPTPl 

05 05 

CE 

18 

IE 


LDK 


»$181E 


0508 

DF 

53 


PRINTS 

STK 


DPTPTR 


0501=1 

CE 

06 

5£ 


LDK 


«HEPD£ 


05 Oil 

BD 

El 

4B 


JSR 


SE14B 

MINIBUS PDPTPl 

051 0 

C6 

05 



LDP 

B 

;s305 


05 IS 

CE 

06 

6D 


LDK 


»HEPD4 


0515 

£0 

03 



ERP 


PRINTS 

. 

0517 

CE 

06* 

73 

PRINT4 

LDK 


«HEPD5 


OSlfl 

BD 


4B 


JSR 


8E14B 

MINIBUS PDPTPl 

05 ID 

CE 

06 

6C 


LDK 


«HEPD3 


05S0 

ED 

El 

4E 

PRINT5 

JSR 


:i;E14B 

MINBUtI PDPTPl 

05S3 

5R 




DEC 

B 



05S4 

£6 

FI 



BNE 


PRINT4 


05S6 

CE 

06 

7E 

LDQP 

LDK 


OHEPD6 

* 

05S9 

BD 

El 

4E 


JSR 


3E14B 

MINIBUS PDPTPl 

05SC 

CE 

00 

4F 


LDK 


JtSPMP 


05SF 

DF 

00 



STK 


KTEMP 


0531 

£0 

£D 



BRP 


LQQPS 


0533 

CE 

06 

tE*'5 

LOOPl 

LDK 


»HEPD7 


0536 

BD 

El 

4B 


JSR 


3E14B 

MINIBUS PDPTPl 

0533 

DE 

00 



LDK 


KTEMP 


053B 

8C 

00 

57 


CPK 


wRDQTPT 


053E 

£6 

IE 



BNE 


LQDP3 


0540 

86 

£0 



LDP 

P 

«3£0 


054S 

ED 

El 

£6 


JSR 


3E126 

QUTCH • 

0545 

EE 

00 



LDK 


OjK 


0547 

E6 

01 



LDP 

E 

1»K 


0543 

FI6 

00 



LDP 

P 

OjK 
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FILTER 5 






054B 

2A 

OB 



BPL 


0S4D 

36 




PSH 

A 

054E 

GE 

06 

SE 


LBX 


0551 

BD 

El 

4B 


JSR 


0554 

32 




PUL 

A 

0555 

BTi 

03 

El 


JSR 


0558 

97 

5F 


LD0P2 

STA 

A 

055A 

D7 

60 



STA 

B 

055C 

20 

06 



BRA 


055E 

EE 

00 


LOOPS 

LBK 


056 0 

EE 

00 


LOOPS 

LBX 


0562 

OF 

5F 



STX 


0564 

BD 

05 

B6 

L0aP4 

JSR 


0567 

BE 

00 



LBX 


0569 

08 




I NX 


056A 

08 




I NX 


0561: 

DF 

00 



STX 


056D 

8G 

00 

59 


CPX 


0570 

26 

Cl 



ENE 


0572 

96 

51 



LDA 

A 

0574 

D6 

52 



LBA 

B 

0576 

CE 

20 



ABB 

B 

0578 

89 

00 



ABC 

A 

057A 

97 

51 



STA 

A 

057G 

D7 

52 



STA 

B 

057E 

96 

53 



LBA 

A 

0580 

D6 

54 



LBA 

B 

0582 

CB 

£0 



ABB 

B 

0584 

89 

00 



ADC 

A 

0586 

97 

53 



STA 

A 

0588 

D7 

54 



STA 

E 

05:5A 

BE 

55 



LBX 


058C 

08 




I NX 


0-53D 

08 




I NX 


058E 

BF 

55 



STX 


0-590 

BE 

57 



LBX 


0592 

08 




I NX 


0593 

08 




I NX 


0594 

BF 

57 



STX 


0596 

BE 

4F 



LBX 


0598 

9C 

59 



CPX 


059A 

27 

OF 



BEQ 


059C 

96 

4F 



LBA 

A 

059E 

B6 

50 



LBA 

B 

OSAO 

CB 

10 



ABB 

B 

05 A2 

89 

00 



ADC 

A 

05A4 

97 

4F 



STA 

A 

05A6 

B7 

50 



STA 

B 

05A8 

7E 

05. 

26 


JMP 



tiECtMBER 19.. 19?8 
LDOPE 
tJMINUS: 

SE14E MI MI BUG PHPTR1 

NEGATE 

Cl 

C2 

L00P4 

0»X 

0»X 

Cl 

COMV 

XTEMP 


XTEMP 

»SftMPMD 

LI30P1 

RAMPTR 

RAMPTR+1 

«S£0 

»;i:0U 

RAMPTR 

RAMPTR-^l 

DATPTR 

DATPTR+1 

j :$20 

:rlrdo 

BfiTPTR 

BATPTR+1 

REPTR 


REPTR 

RBQTPT 


RBOTPT 

SAMP 

SAMPNB 

PEMD 

SAMP 

SAMP+1 

«:i:i 0 

SAMP 

SAMP+1 

LQQP 
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FILTER 5 





DECEMBER 

19» 1973 

OSflB 

©6 

0© 


PEMD 

LB© © 

0© 


05FlD 

C6 

05 



LIi© B 

ttSOS 


aSflF 

BD 

El 

©6 

PEHDl 

JSR 

$E1£6 

MINIBUS OUTCH 






DEC B 



05E3 

©6 

E© 



BNE 

PENDl 


05B5 

3F 




SWI 







♦- 








♦ 




05B$ 

?F 

00 

5E 

COHV 

CLR 

NFLfiS 


05B9 

7D 

00 

5F 

- 

TST 

Cl 


OSBC 

©fi 

03 



BPL 

CONVl 



7© 

00 

5E 


DEC 

NFL©G 


05C1 

CE 

U. 1 

10 

cor+vi 

LBi< 

«3©71 0 

10.- 000 

05C4 

BD 

05 

El 


JSR 

FACTOR 


CISC? 

CE 

03 

E8 


LDX 

«5:03ES 

1> 000 

05Cft 

ID 

05 

El 


JSR 

FACTOR 


05CD 

CE 

00 

64 


LBX 

«30064 

100 

05D0 

BD 

05 

El 


JSR 

FACTOR 


05B3 

CE 

00 

0© 


LDX 

''3000© 

10 

OSDS 

BD 

05 

El 


JSR 

FftCTQR 


GSD9 

36 

60 



LD© ft 

e© 


05DB 

S© 

30 



□R© R 

«330 


05DD 

BD 

El 

©6 


JSR 

3E1 26 

MIN I BUG DUTCH 

G-SEO 

39 




RTS 







♦- 












OSEl 

DF 

61 


FPCTDR 

STX 

D1 


05E3 

?F 

00 

SB 

. 

CLR 

DIOIT 


0516 

ED 

06 

04 

FftCT 

JSR 

SUE2BV 


0519 

©B 

OC 



BMI 

FACTND 


05EB 

?F 

00 

5E 


CLR 

NFLAG 


05EE 

?e 

00 

SB 

F©CT1 

INC 

DIGIT 


05F1 

DE 

63 



LDX 

DIFFl 


05F3 

IF 

5F 



Six 

Cl 


05F5 

20 

EF 



ER© 

FACT 


05F? 

7D 

©0 

5E 

F©CTri© 

TSt 

NFLAG 


05FFI 

©B 

FS 



BNI 

FfiCTl 


06FC 

96 

5B 



LD© © 

DIGIT 


05FE 

©© 

30 



OR© © 

«330 


06 00 

BD 

El 

©6 


JSR 

3E1S6 


0603 

39 




RTS 







♦ 








♦ S BYTE l;UBTR©CTIdN ROUTINE 









06 04 

D6 

60 


SUB2BY 

LD© B 

C© 


0606 

96 

5F 



LD© © 

Cl 

J 

06 0© 

D’O 

6© 



SUB E 

D© 


060© 

9t 

61 



SBC © 

D1 


06 OC 

D7 

64 



ST© B 

DIFF© 
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FILTER 5 


JjECEM£:ER 19.. IS* 78 


06 OE 

97 

63 

0610 

39 


0611 

OD 


0612 

08 

08 

0614 

08 

08 

0616 

20 


0617 

53 

41 

0619 

4D 

50 

06 1 

4C 

‘45 

061 D 

20 

20 

061 F 

20 

20 

062 i 

21 

20 

0623 

20 

20 

0625 

20 

52 

0627 

41 

4B 

0629 

50 

20 

062B 

20 

20 

062D 

20 

21 

062F 

20 


0630 

04 


0631 

20 


0632 

20 

20 


0634 5£ 41 
0636 4Il 50 
063© 20 20 
0638 20 
G63B 04 
063C 20 
0630 20 20 
063F 4E 4F 
0641 49 53 
0643 45 20 

0645 20 

0646 04 

0647 52 
064© 41 40 

0648 50 2B 
064C 4E 4F 
064E 49 53 

0650 45 

0651 04 

0652 20 

0653 21 20 
0655 20 20 
0657 20 20 
0659 52 45 
06 5B 20 20 


ST8 8 DIFFl 
RTS 


♦ 

HE8B 1 FCB 3 OB j 3 08 .» 3 08 .. 3 08 .* 3 08 


FCC S8MPLE • R8MP ! / 


FGB 304 

RHE8D FCC / R8MP / 


FCB 304 

NHE8B FCC / NOISE 


FCB 304 

RNHE8B FCC .'‘R8MP-»^NaiSE/ 


FCB 304 

HE8B2 FCC / ! RE ! RBOT^' 



FILTER 5 


DECEMBER 19-. 1978 


065D 

£0 

£ J 




065F 

£0 

£1 




0S61 

£0 

£0 




0663 

£0 

£0 




0665 

5£ 

44 




0667 

4F 

34 




0669 

01 


HEfllS 

FCB 

SOD? OR? $04 

066R 

OR 

04 




066C 

£D 


HEftD3 

FCC 


066D 

£1 


HER14 

Fce 

\ 

1 

1 

1 

\ 

066E 

£1 

£1 




0670 

£1 

£1 




0678 

£1 

£1 




0674 

£D 

£1 




0676 

£1 

£1 




0678 

04 



FCB 

©04 

0679 

£1 


MER15 

FCE 

©£1>©04 

067R 

04 





067E 

00 

01 

HER16 

FIB 

© 01 j © OR !» ©£ 0 j ©£ 0 !. © 0 4 

067D 

00 

OR 

n 




1 P 

0681 

uu 

00 

e U 
20 


- 


0683 

00 

04 




0685 

£0 


MER17 

FCC 

/ ! / 

06w6 

£0 

£0 




068© 

£0 

£1 




06SR 

£0 

£0 




06SC 

£0 




• 

06©B 

04 



FCB 

©04 

068E 

0 © 


Minus 

FCB 

©OS!.©£D»©04 


06©F 81 04 

EMI 

MO ERROR <:S> DETECTED 


SYMBOL TRBLE: 


R1 

OO'OE 

R£ 

OOOF 

R3 

0010 

Rll 

03E9 

B1 

0011 

B£ 

0O1£ 

B3 

0013 

Cl 

005F 

C£ 

0060 

CMT 

0003 

ClM'v' 

0SB6 

GOMVl 

05C1 

11 

0061 

D£ 

O06£ 

IRTPTR 

0053 

DIFFI 

0063 

1IFF£ 

0064 

IISIT 

0O5B 

EMI 

Oc*c*F 

EMIFL© 

0067 

FRCT 

05E6 

FRCTl 

05EE 

FRCTM® 

05F7 

FRGT0R 

OSEl 

FLT 

0200 

FLTl 

0256 

FLT 11 

U2S0 

FLT£ 

0£95 

FLT3 

0£9R 

FLT4 

O£E0 

FLTR 

0£41 

NERll 

061 1 

WER1£ 

065£ 

HERD3 

066C 

HER14 

0661 

WER15 

0679 



FILTER 5 


IiECEMUER 19 f 197S 


HERDS 

UST*B 

HERB? 

0335 

MER03 

0339 

INITl 

04 IF 

I NITS 

0431 

INPUTl 

0223 

INPUTS 

OtSE 

INPUTS 

0231 

LOOP 

0526 

LOOPl 

0533 

LQDP2 

0553 

LOOPS 

055E 

L00P4 

05S4 

L0QP5 

0530 

MINUS 

06SE 

MLTCDl 

0004 

MLTC02 

0005 

MLTPRl 

0006 

MLTPR2 

000? 

MLTPR3 

0 0 03 

MULT 

0347 

MULTI 

035F 

MULT2 

037? 

MULT3 

0370 

MULT31 

0330 

MULT4 

0397 

MULT5 

0303 

MULTND 

0309 

NEGl 

.03E? 

NEG2 

03E3 

NEG3 

O30C 

NEG4 

03E2 

NEGATE 

03E1 

NFLRG 

0O5E 

NHER0 

0333 

OUT 

03 IB 

OUTl 

03S9 

PEN0 

osr.B 

PEN01 

05RF 

PR INTI 

04F1 

PRINTS 

04FF 

PRINT3 

0503 

PRINT4 

051? 

PRINTS 

0520 

PRINTR 

04C7 

pRnii 

0009 

PRD02 

OOOR 

PRO03 

OOOB 

PRDD4 

OOOC 

PRQBS 

0000 

RRMPTR 

0051 

RCNT 

0044 

RCRCtl 

0043 

RGRCT2 

0047 

RDI 

0OS5 

RBOTl 

0043 

R0OT2 

0049 

RDOTPT 

0057 

REl 

004R 

RE2 

0O4B 

REGS 

0020 

REGSN0 

0040 

REPTR 

0055 

RMERD 

033 1 

RNHEP0 

0347 

R0UN0 

0050 

RDUN01 

039B 

ROUN0C’ 

03E2 

RTEMP 

001 E 

SRMP 

O04F 

SRMPN© 

0059 

SGN 

0002 

SHFTR 

04 OF 

SHFTR 1 

0414 

SKPFLG 

O'045 

SUB 

03FC 

SUB’EBY 

0304 

SUMl 

0014 

SUMS 

0015 

SUMS 

0016 

TSUMl 

O04G 

TSUMS 

0040 

TSUM3 

004E 

TYPE 

OOSC 


XTEMF 0000 



APPENDIX D 


The program to determine variance reduction ratio from the 
2 transfer function is presented on the following page. 
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PROGRAM TO DETERMINE VARIANCE REDUCTION RATIO FROM Z TRANSFER FUNCTION 


1 0 

IiIf’IN C5 

15 

S=0 

£0 

PF:INT‘'I 

£3 

INPUT ■ 

4 0 

PRINT”! 

50 

INPUT 1 

&0 

FDR J= 

70 

PRINT"! 

3 0 

INPUT ! 

30 

MEKT J 

1 0 0 

PRINT"! 

no 

INPUT : 

ISO 

FDR J= 

1 3 0 

PRINT": 

140 

INPUT . 

150 

NEXT J 

155 

PRINT 

160 

IF R>=, 

170 

R=5‘ 

1 S 0 

FDR J=^l 

1 3 0 

N < J> =0 

£ 0 0 

NEXT J 

SI 0 

FDR L = 

£20 

M=N < 1 . 

£30 

FDR J=: 

£40 

N CJ;:< =H 

£5 0 

NEXT J 

£6 0 

S=3+M'f- 

£70 

PRINT" 

£3 0 

FOR 

£'3 0 

N CJ;:' =N 

3 0 0 

NEXT J 

3 1 0 

N <Pi-!-i;5 

330 

NEXT L 

36 0 

EfiJJ 


R>=B THEN £10 


T-i i \ 




APPENDIX E. 


A method for deriving a discrete system from a continuous system 

is to replace the integrators in the continuous system with numerical 
9 

integrators. The resulting discrete system will be different for each 
method of numerical integration used. The simplest method is Euler's 
approximation where the time derivative is approximated by 


|_f(„T) . m) - f»n-l)T) . 


(37) 


Taking the Laplace operator of the left side of (37) and equating it to 
the z operator of the right side of (37) yields the mapping from the 
s-plane to the z-plane 


l-z“ 

t ' • 


(33) 


The continuous integrator can then be replaced by 


i - t 

sr. " I-z^^ 



( 39 ) 



203 


Although with this method a discrete system can be obtained for 
any continuous system, the discrete system may not be physically realizable 
on a digital computer. Due to the finite time required for computations, 
every feedback loop must contain a delay element. 


